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SUMMARY 


This  report  develops  models  for  determining  collision-limiting,  or  “delta” 
factors  for  certain  electronic  states  of  iVj,  N^,  Oj,iVO  and  NO^  molecules.  These 
factors,  when  multiplied  by  the  Boltzmann  populations  of  the  respective  states,  yield 
the  steady-state  populations  as  they  exist  under  the  influence  of  electron  and  heavy- 
pajticle  collisions  and  a  radiation  field.  With  one  exception,  only  electric-dipole- 
allowed  transitions  are  included.  The  resulting  radiation  lies  mainly  in  the  visible 
and  ultraviolet  (UV)  regions.  Since  the  volume-emission  rate  is  proportional  to  the 
population  of  the  emitting  state,  the  delta  factors  can  also  be  used  to  convert,  directly, 
thermal  radiation,  calculated  under  conditions  of  local  thermodynamic  equilibrium 
(LTE)  at  the  neutral  gas  temperature,  to  that  existing  under  conditions  of  non-local 
thermodynamic  equilibrium  (NLTE). 

For  each  molecule,  and  for  each  state  of  each  molecule  considered,  a  steady- 
state  population  is  determined  from  a  balance  between  its  rate  of  population  and 
depopulation  by  collisional  and  radiative  processes.  Separate  electron  and  neutral-gas 
temperatures  are  provided  for.  Radiative  and  collisional  coupling  among  the  states 
of  a  given  species,  where  significant,  necessitates  the  solution  of  simultaneous  linear 
equations  for  the  state  populations  involved.  For  radiative  absorption,  only  sunlight 
is  included  in  our  models. 

The  number  and  type  of  states  for  which  delta  factors  are  determined,  and 
the  radiative  band  system'^  involved,  are  as  follows.  For  iVj,  there  are  11  electronic 
states  including  the  ground.  The  band  systems  are  the  Vegard-Kaplan,  First  Positive, 
Second  Positive,  Birge-Hopfield  1,  Birge-Hopfield  2  and  five  different  Rydberg  series. 
For  Nj,  there  are  4  electronic  states  including  the  ground.  The  bands  are  the  Meinel, 
First  Negative  and  Second  Negative  Systems.  For  0%,  besides  the  ground  state,  we 
include  21  vibrational  states  of  a  single  electronic  state.  Transitions  from  these  states 
to  the  ground-vibrational  states  give  rise  to  bands  of  the  Schumann-Runge  System. 
For  NO,  there  au-e  7  electronic  states  including  the  ground.  The  band  systems  are  the 
gamma,  beta,  delta.  Heath,  epsilon.  Feast  1,  Feast  2,  beta  prime,  gamma  prime,  Feast- 
Heath  and  —  C*n.  Finally,  for  NO'*',  deltas  are  determined  for  the  ground  state 
and  for  10  vibrational  states  of  a  single  electronically-excited  state.  Transitions  from 
these  states  to  the  ground-  vibrational  states  give  rise  to  bands  of  the  Baer-Miescher 
System  in  the  vacuum  UV. 


The  models  have  been  computerized  with  a  main  diver  program  DELTAVU 
and  a  separate  subroutine  for  each  of  the  5  molecules  involved.  Sample  calculations 
are  shown  giving  delta  factors  computed  for  three  different  nuclear-burst  conditions 
as  determined  by  the  NORSE  code. 

Some  recommendations  are  made  regarding  future  improvements  and  pos¬ 
sible  expansion  of  the  models  to  include  additional  states  and  species. 
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CONVERSION  TABLE 


Conversion  factors  for  U.S.  Customary  to  metric  (SI)  units  of  measurement 

MULTIPLY  - -  BY  - -  TO  GET 

TO  GET  - -  BY  - -  DIVIDE 


angstrom 

1.000000  X  E  -10 

meters  (m) 

atmosphere  (normal) 

1.01325  X  E  -1-2 

kilo  pascal  (kPa) 

bar 

1.000000  X  E  -f-2 

kilo  pascal  (kPa) 

barn 

1.000000  X  E  -28 

meter^  (m’) 

British  thermal  unit  (thermochemical) 

1.054350  X  E  -t-3 

Joule  (J) 

calorie  (thermochemical) 

4.184000 

joule  (J) 

cal  (thermochemical)  /  cm^ 

4.184000  X  E-2 

mega  joule/m'  (MJ/m^) 

curie 

3.700000  X  E  -^l 

•giga  becquerei  ’’•Bq) 

degree  (angle) 

1.745329  X  E-2 

radian  (rad) 

degree  Farenheit 

‘jf  =  (‘f  +  459.67)/1.8 

degree  kelvin  (K) 

electron  volt 

1.60219  X  E-19 

joule  (J) 

erg 

1.000000  X  E  -7 

joule  (J) 

erg/second 

1.000000  X  E-7 

watt  (W) 

foot 

3.048000  X  E -1 

meter  (m) 

foot-pound-force 

1.355818 

joule  (J) 

gallon  (U.S.  liquid) 

3.785412  X  E  -3 

meter*  (m^) 

inch 

2.540000  X  E  -2 

meter  (m) 

jerk 

1.000000  X  E  -t-9 

joule  (J) 

joule/kilogram  (J/kg)  (radiation  dose  absorbed) 

1.000000 

Gray  (Gy) 

kilotons 

4.183 

terajoules 

kip  (1000  Ibf) 

4.448222  x  E  -f3 

newton  (.V) 

kip/inch^  (ksi) 

6.894757  x  E  4-3 

kilo  pascal  (kPa) 

ktap 

1.000000  X  E  -f-2 

newton-second/ m’  (N-s/m’) 

micron 

1.000000  X  E  -6 

meter  (m) 

mil 

2.540000  X  E  -5 

meter  (m) 

mile  (international) 

1.609344  X  E  4-3 

meter  (m) 

ounce 

2.834952  x  E  -2 

kilogram  (kg) 

pound-force  (lbs  avoirdupois) 

4.448222 

newton  (N) 

pound-force  inch 

1.129848  X  E-1 

newton-meter  (N  m) 

pound- fcrce/inch 

1.751268  X  E  4-2 

newton/meter  (N/m) 

pound-force/foot’ 

4.7S8026  X  E  -2 

kilo  pascal  (Kra; 

pound-force/inch^  (psi) 

6.894757 

kilo  pascal  (kPa) 

pound-mass  (Ibm  avoirdupois) 

4.535924  x  E  -1 

kilogram  (kg) 

pound-mass-foot^  (moment  of  inertia) 

4.214U11  X  E -2 

kilogram-meter’  (kgm’) 

pound-mass/foot* 

1.601846  X  E  4-1 

kilogram/ meter’  (kg/m’) 

rad  (radiation  dose  absorbed) 

1.000000  X  E  -2 

••Gray  (Gy) 

roentgen 

2.579760  X  E  -4 

coulomb/kilogram  (C/kg) 

shake 

1.000000  X  E  -8 

second  (s) 

slug 

1.459390  X  E  4-1 

kilogram  (kg) 

torr  (mm  Hg,  0"  C) 

1.333220  X  E-1 

kilo  pascal  (kPa) 

*The  becquerel  (Bq)  ia  the  SI  unit  of  radioactivity;  1  Bq  =  1  event/s. 
**The  Gray  (Gy)  i»  the  SI  unit  of  absorbed  radiation. 


VI 


TABLE  OF  CONTENTS 


Section  Page 

SUMMARY  .  iii 

PREFACE  .  V 

CONVERSION  TABLE .  vi 

LIST  OF  TABLES .  ix 

1  INTRODUCTION .  1-1 

2  DELTA  MODEL  FOR  Nj .  2-1 

STATES  AND  BAND  SYSTEMS  INCLUDED .  2-1 

EXCITATION/DEEXCITATION  MECHANISMS .  2-2 

DELTA-FACTOR  EQUATIONS  FOR  N, .  2-6 

Singlet  States .  2-7 

Triplet  States .  2-9 

Ground-State  Delta  Factor,  .  2-10 

Electronic  Partition  Function  for  Nj .  2-11 

3  DELTA  MODEL  FOR  NJ^  .  3-1 

STATES  AND  BAND  SYSTEMS  INCLUDED .  3-1 

EXCITATION/DEEXCITATION  MECHANISMS .  3-1 

DELTA-FACTOR  EQUATIONS  FOR  N^ .  3-5 

4  DELTA  MODEL  FOR  Oj .  4-1 

STATES  AND  BAND  SYSTEMS  INCLUDED .  4-1 

EXCITATION/DEEXCITATION  MECHANISMS .  4-1 

Electron  Excitation/Deexcitation  Rate  Coefficients .  4-2 

vii 


TABLE  OF  CONTENTS  (Concluded) 

Section  Page 

Einstein  A  Coefficients  and  pB  Factors  .  4-4 

DELTA-FACTOR  EQUATIONS  FOR  Oj .  4-4 

5  DELTA  MODEL  FOR  NO .  5-1 

STATES  AND  BAND  SYSTEMS  INCLUDED .  5-1 

EXCITATION/DEEXCITATION  MECHANISMS .  5-1 

Electron-Impact  Excitation/Deexcitation .  5-5 

Heavy-Particle-Impact  Excitation/Deexcitation .  5-6 

Radiation  Parameters .  5-7 

DELTA-FACTOR  EQUATIONS  FOR  NO .  5-7 

States  B*n  and  5-8 

States  and  B*E+ .  5-9 

6  DELTA  MODEL  FOR  NO+  .  6-1 

STATES  AND  BAND  SYSTEMS  INCLUDED .  6-1 

EXCITATION/DEEXCITATION  MECHANISMS .  6-1 

Electron  Excitation/ Deexcitation  Rate  Coefficients .  6-2 

Einstein  A  Coefficients  and  pB  Factors .  6-3 

DELTA-FACTOR  EQUATIONS  FOR  NO^ .  6-4 

7  SOME  NUMERICAL  EXAMPLES .  7-1 

8  COMMENTS  AND  RECOMMENDATIONS .  8-1 

9  REFERENCES .  9-1 

Appendix 

FORTRAN  LISTING  OF  PROGRAM  DELTAVU  AND 
ASSOCIATED  SUBROUTINES  .  A-1 


viii 


LIST  OF  TABLES 

Table  Page 

2.1  States  of  Nj  included  .  2-1 

2.2  Band  systems  of  Nj  included .  2-2 

2.3  Excitation/deexcitation  mechanisms  and  assumed  rates  for 

states  of  N  j  .  2-4 

2.4  Delta-factor  notation  for  states  of  Nj  .  2-7 

2.5  Definition  of  parameters  in  Equation  2.1.* .  2-10 

3.1  States  of  Nj  included .  3-1 

3.2  Band  systems  of  NJ”  included .  3-2 

3.3  Excitation/deexcitation  mechanisms  and  assumed  rates 

for  NJ’ .  3-3 

3.4  Delta-factor  notation  for  states  of  NJ" .  3-5 

4.1  State  and  band-system  parameters  for  Oj .  4-1 

5.1  States  of  NO  included .  5-1 

5.2  Band  systems  of  NO  included .  5-2 

5.3  Excitation/deexcitation  mechanisms  and  assumed  rates  for  states  of 

NO .  5-3 

5.4  Average  excitation  energy  auid  electronic-absorption  oscillator  strengths 

for  NO  transitions .  5-6 

5.5  Delta-factor  notation  for  states  of  NO .  5-7 

6.1  States  and  band-system  parameters  for  NO''" .  6-1 

7.1  Delta  factors  at  1  sec  in  a  cell  at  200-km  altitude,  horizontally  dis¬ 
placed  140  km  from  a  MT-class  burst  at  200-km  altitude  ....  7-4 

7.2  Delta  factors  at  1  sec  in  a  cell  at  140-km  altitude,  directly  above  a 

MT-class  burst  at  90-km  altitude .  7-5 

7.3  Delta  factors  at  1  sec  in  a  cell  at  140-km  altitude,  horizontally  dis¬ 
placed  50  km  from  a  MT-class  burst  at  600-km  altitude .  7-6 


IX 


SECTION  1 


INTRODUCTION 


This  report  is  the  visible/ultraviolet  (UV)  counterpart  to  References  1.1  and 
1.2  that  describe  a  model  for  computing  the  non-local  thermodynamic  equilibrixun 
(NLTE)  radiation  from  selected  IR-active  molecules.  In  those  reports,  we  dealt  with 
the  determination  of  non-Boltzmann  populations  of  vibrational  states  that  give  rise  to 
the  IR  radiation.  In  this  report,  we  are  interested  in  determining  the  non-Boltzmann 
populations  mainly  of  electronic  states  of  the  air  species  iVj,  Nj  ,  Oj,  NO  and  NO'*' 
that  give  rise  to  visible/UV  radiation.  The  object  is  the  construction  of  NLTE  emission 
models  for  incorporation  in  the  NORSE  code. 

The  approach  taken  here,  as  in  References  1.1  and  1.2,  is  one  of  providing 
a  collision-limiting,  or  “delta”  factor  for  each  radiating  state  considered.  This  factor 
is  defined  as  one  which,  when  multiplied  by  the  Boltzmann  population  of  the  state 
(or  by  the  LTE  emission  rate  from  it),  will  yield  the  corresponding  non-Boltzmann 
population  (or  the  NLTE  emission  rate).  The  method  entails  the  calculation  of  the 
steady-state  concentration  of  each  state  by  equating  its  population  and  depopulation 
rates  from  collisional  and  radiative  processes. 

The  electronic  states  chosen  for  inclusion  in  our  model  are  those  in  the  Air 
Force  Weapons  Laboratory  (AFWL)  opacity  model  (Reference  1.3).  Generally,  only 
permitted  transitions  (i.e.  electric-dipole  allowed)  involving  these  states  are  included. 
This  is  consistent  with  the  NORSE  thermal-emission  model  that  uses  absorption  cross 
sections  for  these  same  transitions  (bands)  in  determining  the  LTE  volume  emission 
rates  for  the  visible/UV. 

Except  for  the  cases  of  Oj  and  NO"*",  for  which  only  one  electronically- 
excited  state  is  significant,  a  single  delta  factor  is  modeled  for  each  electronic  state 
considered.  That  is,  with  the  above  exceptions,  the  same  delta  factor  is  assumed  to 
apply  to  each  vibrational  state  of  a  given  electronic  state.  Thus,  the  same  delta  factor 
is  used  to  convert  the  LTE  emission  from  all  bands  of  a  given  band  system  to  NLTE 
emission.  This  is  done  for  three  reasons.  First,  NORSE  cannot  presently  identify 
contributions  from  each  band  of  a  given  band  system  to  the  LTE  emission  in  a  spec¬ 
ified  wavenumber  interval.  Therefore,  it  is  unable  to  utilize  the  increased  resolution 
provided  by  the  availability  of  separate  deltas  for  each  band  of  a  given  band  system. 
Second,  insufficient  data  exist  on  excitation/ deexcitation  rates  for  the  individual  vi¬ 
brational  states  to  permit  development  of  a  reliable  delta-factor  model  for  a  multitude 
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of  vibrational  states  of  each  electronic  state.  Third,  the  development  of  such  a  model 
would  be  quite  unwieldy,  even  if  adequate  data  were  available.  Consequently,  the 
models  developed  in  this  report  determine,  in  some  sense,  average  delta  factors  for 
entire  band  systems  and,  as  such,  they  provide  what  can  be  considered  first-order 
corrections  to  the  LTE  emission  from  each  band  of  a  given  system. 

Coupling  between  electronic  states  of  different  molecules  is  assumed  to  be 
negligibly  small  and  so  a  separate  delta  model  can  be  developed  for  each  molecule, 
independent  of  that  for  each  of  the  others.  Sections  2  through  6  below  describe 
the  delta  models  for  iVj,  N^,  Oj,  NO  and  iVO"*",  respectively.  Section  7  shows  sam¬ 
ple  delta-factor  results  and  Section  8  presents  some  comments  and  recommendations. 
The  models  have  been  computerized  with  a  main  driver  program  DELTAVU  and  sepa¬ 
rate  subroutines  DELN2,  DELN2P,  DEL02,  DELNO  and  DELNOP  for  the  molecules 
jVj,  ,  O2,  NO  and  NO'*' ,  respectively.  The  computer  listings  are  given  in  the  Ap¬ 
pendix. 
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SECTION  2 


DELTA  MODEL  FOR  Nj 


STATES  AND  BAND  SYSTEMS  INCLUDED 


Table  2.1  shows  the  electronic  states  of  Ni  included  in  our  Delta  model. 
Also  shown  are  the  energy  (Ref. 2.1)  and  statistical  weight  of  each  state.  The  energy 

Table  2.1.  States  of  N]  included. 


Statistical  Weight 

g 

o.o/o.o 

1 

49,754.78/6.18 

3 

J5^n, 

59,306.81/7.37 

6 

c»n„ 

88,977.84/11.05 

6 

n„ 

100,817.5/12.53 

2 

103,670.8/12.86 

1 

cz  ^n„(c  ^n„) 

104,138.2/12.92 

2 

c'.  ■e;(c'  'e;) 

104,323.3/12.95 

1 

oz  ‘n„(o  ‘n„) 

105,710.4/13.12 

2 

C4  ‘n„(e  ‘n„) 

115,635.9/14.34 

2 

115,876/14.37 

1 

values  shown  are  for  the  ground  vibrational  states.  The  statistical  weight,  g,  for  E 
states  is  2S  +  1  while  for  other  states  (11,  A,  •••)  it  is  2(2S+1)  (Ref.  2.2)  where 
5  is  the  total  spin  quantum  number.  The  first  six  entries  in  Table  2.1  are  valence 
states;  the  last  five  entries  are  Rydberg  states.  For  the  Rydberg  states,  quantities 
in  parentheses  indicate  the  older  state  notation;  numerical  subscripts  refer  to  the 
principal  quantum  number,  n.  The  state  os  is  the  first  member  (n  =3)  of  the 
Rydberg  series  o„^n(n  =  3,4,5,---)  converging  on  iV^(A*n„),  while  the  other  four 
Rydbergs  are  members  of  series  converging  on  iVj'^(X’Ej). 


Table  2.2  shows  the  band  systems  of  iVj  included  In  our  model.  All  are 
permitted  transitions  except  for  the  Vegard-Kaplan  System  that  violates  the  AS  = 
0  rule  for  electric-dipole  radiation.  Also  shown  in  Table  2.2  are  the  (0-0)  band-origin 
wavelengths  and  wavenumbers  of  the  systems  as  well  as  the  approximate  range  of 
wavelengths  over  which  each  band  system  radiates  (Ref.  2.1). 

Table  2.2.  Band  systems  of  Nj  included. 


Baind  System 

Aqo 

(A) 

I'OO 

(cm-‘) 

Wavelength  Region 

(A) 

(Vegard-Kaplan) 

2,009.9 

49,754.8 

5325-1250 

(First  Positive) 

10,469.0 

9,552.0 

25,310-4780 

— B^Ug  (Second  Positive) 

3,370.3 

29,671.0 

5460-2680 

6^n„— (Birge-Hopheld  1) 

991.9 

100,817.5 

1110-850 

6'  ‘E+— X‘E;  (Birge-Hopfield  2) 

964.6 

103,670.8 

1290-820 

C3  (Worley-Jenkins  Ryd.) 

960.3 

104,138.2 

960-922 

c;  (Carroll- Yoshino  Ryd.) 

958.6 

104,323.3 

960-860 

03  Ar‘E;  (Worley’s  Third  Ryd.) 

946.0 

105,710.4 

950-880 

c^  — X'E^  (Worley-Jenkins  Ryd.) 

864.8 

115,635.9 

-865 

Cj  ‘Ey — X‘E^  (Carroll- Yoshino  Ryd.) 

863.0 

115,876 

-863 

EXCITATION /DEEXCITATION  MECHANISMS 


The  processes  that  we  include  for  excitation  auid  deexcitation  of  the  elec¬ 
tronic  states  are  electron  and  heavy-particle  collisions,  spontaneous  radiation  and 
sunlight  absorption.  The  information  that  we  have  on  the  collisional  processes  is  lim¬ 
ited.  For  the  singlet  states  considered  here,  we  have  found  no  data  on  excitation  or 
quenching  by  heavy  particles  and  so  we  include  these  processes  only  for  electrons. 
However,  this  should  not  be  a  serious  limitation  since,  for  a  given  temperature,  the 
electron  processes  usually  proceed  at  much  faster  rates  than  do  those  for  heavy  par¬ 
ticles  (Ref.  2.3).  For  the  triplet  states,  data  on  heavy-particle  collisions  is  limited  to 
room  temperature  or  to  narrow  temperature-range  measurements. 
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Many  of  the  high-lying  singlet  states  (valence  and  Rydberg)  are  strongly  pre¬ 
dissociated  and  their  radiation  is  thus  diminished  (Refs.  2.1,  2.4).  Although  sufficient 
information  is  available  on  the  predissociation  rates  to  make  approximate  allowance 
for  them,  the  reverse  process,  preassociation,  that  will  tend  to  partially  compensate 
(Ref.  2.5),  cannot  be  so  easily  allowed  for.  For  this  reason  we  ignore  these  processes 
in  our  model.  As  a  consequence,  the  delta  factors  predicted  by  our  model  for  these 
states,  and  thus  the  NLTE  emission  rates  from  them,  tend  to  be  upper  limits  to  the 
true  values. 

Table  2.3  shows  the  specific  reactions  or  processes  included  in  the  iVj  model 
together  with  suggested  values  for  the  rate  coefficients,  the  Einstein  A  coefficients  and 
pB  factors  ^  for  sunlight.  The  forward  (deexcitation)  rate  coefficients  for  the  first  19 
reactions  in  Table  2.3  are  expressed  in  the  form 


k  =  (cm*  sec~^)  (2.1) 

where  the  parameters  a,  6  and  c  are  given  in  Columns  2,  3  and  4,  respectively,  while 
the  reverse  (excitation)  rate  coefficients  are  expressed  in  the  form 

k' =  (cm*sec-‘)  (2.2) 

where  a'  and  c'  are  given  in  Columns  5  and  6,  respectively.  For  collisions  involving 
heavy  particles,  T  is  the  heavy-particle  temperature  (°K)  while,  for  collisions  involving 
electrons,  it  refers  to  the  electron  temperature  T*.  Also  shown  in  Table  2.3  are  the 
temperatures,  or  temperature  ranges,  for  which  the  rate-coefficient  expressions  apply 
and  the  supporting  references. 

With  respect  to  Reactions  2-11  through  2-19,  we  should  point  out  that  the 
products  on  the  right-hand  sides  au'e  generally  uncertain,  although  there  is  mounting 
evidence  that  the  product  state  of  Wj  is  usually  not  the  ground  state,  espe¬ 

cially  where  quenching  of  the  B  and  C  states  is  concerned.  Indeed,  the  reactions  can 
be  chemical  in  nature  with  the  product  species  being  other  than  N2.  Nevertheless, 
for  the  sake  of  simplicity  in  allowing  for  the  reverse  reactions,  we  assume  that  the 
electronically-excited  Nj  molecules  are  quenched  to  the  ground  state. 

'■p  is  the  energy  density  of  the  radiation  field  at  the  band  frequency  and  B  is  the  Einstein  B 
coefficient.  The  product  has  dimensions  8ec~^. 
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Table  2.3.  Excitation/deexcitation  mcchaniBrns  and  assumed  rates  for 
states  of  N]. 
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Tabic  2.3.  Excitation/decxcitation  mechanums  and  assumed  rates  fur 
states  of  TVj  (Continued). 


The  rate  coefficients  for  excitation  by  electron  impact  (Reactions  2-1  through 
2-10  in  Table  2.3)  were  determined  from  the  measured  cross  sections,  <7{E),  of  Cartwright 
et  al.  (Ref.  2.6)  using  the  relation  (Ref.  2.7) 

k'  =  6.69  X  10-®rf /“  <t{E)E  e-^f^’dE  .  (2.3) 

JEo 

Here,  the  electron  energy,  E,  the  threshold  energy  for  excitation..  Eq,  and  the  electron 
temperature,  T*,  are  all  in  eV  units  and  tT{E)  is  in  A’.  We  performed  the  integrals 
in  Equation  2.3  numerically  for  parametric  values  of  T*  ranging  from  0.1  to  3.0  eV. 
The  results  were  fitted  to  the  form  2.1  with  T*  in  “K  and  they  were  then  reexpressed 
in  the  form  2.2  using  the  relationship  between  the  excitation  and  deexcitation  rate 
coefficients: 


jfe'  = 

Qm 


(2.4) 


where  En  >  E^  and  fcj  is  Boltzmann’s  constant. 

Reactions  2-20  through  2-29  in  Table  2.3  are  the  radiation  processes  included 
in  the  model.  For  these,  Columns  2  and  3  give  the  Einstein  A  coefficients  (sec~^)  and 
the  pB  factors  (sec~^)  for  sunlight  conditions,  respectively.  The  latter  are  based  on 
the  relation  (Ref.  1.1) 


pB  =  5.4  X  10-«  (e2.39ixio-«p„.,  _  ,2.5) 

9m  ^  ’ 

applicable  for  a  6000  "K  sun  subtending  an  amgle  of  6.8  x  10“®  ster  at  the  earth. 
The  Einstein  coefficients  in  Table  2.3  are  “generic”  in  the  sense  that  no  one  value 
corresponds,  necessarily,  to  precisely  the  value  for  a  specific  band  of  a  given  band 
system  but,  rather,  is  an  average  of  the  values  in  the  referenced  data  for  various 
bands  of  the  band  system. 


DELTA-FACTOR  EQUATIONS  FOR  Nj 

As  shown  in  Table  2.1,  there  are  eleven  states  to  determine  delta  factors  for. 
Table  2.4  shows  the  notation  used  for  these  factors. 
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Table  2.4.  Delta-factor  notation  for  states  of  N2. 


State 

Delta 

State 

Delta 

State 

Delta 

dj'* 

Os  ‘n,. 

wm 

df* 

cs  ‘Hu 

•s’e: 

dl?^ 

From  the  set  of  29  reactions  in  Table  2.3,  it  is  clear  that  in  our  model  the 
seven  excited  singlet  states  e^e  coupled  only  to  the  ground  state,  whereas  the  three 
triplet  states  are  coupled  to  each  other  as  well  as  to  the  ground  state.  Determination 
of  the  NLTE  populations  and  delta  factors  for  the  singlet  states,  therefore,  is  simpler 
than  that  for  the  triplet  states. 


Singlet  States 

As  an  example,  we  now  illustrate  the  procedure  for  determining  the  delta 
factor  for  the  6‘nu  state.  Determination  of  the  delta  factors  for  the  other  singlet 
states  follows  in  a  similar  manner. 

By  reference  to  the  reaction  set  in  Table  2.3,  the  time  rate  of  change  of  the 
concentration  of  state  6‘nu  can  be  written  as 


=  {k[N,  +  -  [kiN,  +  A2o){N,{b^U^)\  .  (2.6) 

The  steady-state  solution  is 


where  iV«  is  the  electron  density  and  subscripts  on  k,  k',  A  and  pB  correspond  to  the 
reaction  numbering  in  Table  2.3.  By  definition,  the  delta  factor  for  state  is  given 
by  the  ratio 
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(2.8) 


where  the  denominator  is  the  Boltzmann  population  of  the  state.  But  the  Boltzmann 
population  for  state  n  of  an  arbitrary  species  V  can  be  expressed  in  the  form 

[y{n)U  =  (2.9) 

where  [V]  is  the  species  concentration,  En  is  the  state  energy  and  Qe(T'),  the  electronic 
partition  function,  is  defined  as 


Qc{T)  = 


n=0 


(2.10) 


Alternatively,  in  terms  of  the  ground-state  concentration  [3^(0)j,  Equation  2.9  can  be 
written 


|V'(n)W  =  -  =  ^|y(0)U.-‘  (2.11) 

9o  go 

where  i/„  is  the  energy  (cm“^)  of  state  n.  Therefore,  given  the  energy  and  statistical 
weight  in  Table  2.1  for  state  6^n„,  it  follows  that 


[iV,(6‘nu)U  = 

Equations  2.7,  2.8  and  2.12  now  lead  to  the  result 

Here,  is  the  delta  factor  for  the  ground  state: 

=  [;v,(A:‘E,^)i„/(iV:(x‘E;)i.,  . 

Determination  of  will  be  addressed  below. 


(2.12) 


(2.13) 


(2.14) 
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In  a  similar  fashion,  we  arrive  at  the  following  expressions  for  the  delta 
factors  for  the  remaining  singlet  states: 


(  k'jNt  +  pBji  ^  i49iao.j/r  jyvj 

1  /  k^Ng  4-  pBn  \  i49a3j.7/r^j 

2  kiNg  +  J  ^ 

I  k\Ng  +  pBis\  150099.0/T 

[  k,Ng  +  y  ° 

1  /  k'^Ng  +  (>Bn\  I5j094.7/r^a 

2  I,  k,Ng  +  A2J  ^ 

1  I  k'^Ng  +  pB2s\  ^166375.4/r^a 

2  \  k^Ng  +  Ais  } 


“10 


(k'^Ng+pB^A 

V  ^7^e  +  -^18  / 


gl667a0.9/T^^j 


(2.15) 

(2.16) 


(2.17) 


(2.18) 

(2.19) 

(2.20) 


Triplet  States 

As  indicated  above,  the  three  triplet  states  are  coupled  to  eau:h  other  and 
to  the  ground  state  by  Reactions  2-8  through  2-19  and  2-27  through  2-29  in  Table 
2.3.  In  steady  state,  the  delta  factor  for  each  triplet  can  then  be  written  as  a  linear 
combination  of  the  deltas  for  the  other  two  triplets  and  for  the  ground  state.  In 
particular,  we  find  that  the  delta  factors  df*,  and  are  related  to  d^’,  the 
ground-state  delta  factor,  by  meams  of  the  equations 


ajdJ'*  -  aad^’  +  O.d^^  =  a^d^^ 

P2d^^  -  03d^^  +  P,d^^  = -0,d^^  (2.21) 

O.df^  +  isd^^  -  l*d^^  = -lid^^  , 


where  the  coefficients  are  defined  in  Table  2.5. 
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Table  2.5.  Definition  of  parameters  in  Equation  2.1.*. 


0(1  —  ^16,17,18,19(^/1  + 

—  ^is.w.isl-^l  +  ^9-^e  +  -^28 

aj  =  ^16,1T.18,19(-^)  “»■  klQ^t  +  -^29 

^4  =  : 

as  =  224s8er-i”«^/^ 

'/i  =  ikiU2[A/]  +  fc8iV,e‘”‘”°^‘r-2ir)  . 

/3i  =  A:i3,i4,i6lA/]  + 

7s  » 

3,  = 

l4  =  ^n,12[A/]  +  +  24j7 

1  *  The  notation  ^18,17,18, i9[A/j  is  short  for  ^lelOj)  -1-  l:i7[iVO]  -t-  A:i8[G]  -1-  l:i9[.V|. 

I  Similarly,  Icis  14  islM]  ^11,12!^]  follow  from  Table  2.3. 

Equations  2.21  are  readily  solved  for  and  d^'‘  to  yield 

<i\‘‘  =  +  ^3^'1x)\d^^  /  D 

=  Ml<  +  M0nA+/34lxM^/D  (2.22) 

d^^  =  [ai/?i73  +  or2(/?3Tfi  +  Pils)  -  az0,'ii]d^^ / D 

where 

D  =  ai{0sni  -  04lz)  -  azPilA  ■  (2.23) 

It  remains  to  determine  the  ground-state  delta  factor,  . 

Ground-State  Delta  Factor, 

In  general,  for  an  arbitrary  species  Y  with  electronic  states  n  =  0, 1,2, - . 

we  can  write 

!V'(o)]  +  iy{i)i  +  iy(2)|  + - =  iy|  .  (2.21) 
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Dividing  both  sides  of  the  Equality  2.24  by  [^^(O)]^,,  and  recognizing  that 


ywi  ^  |ywi  i^wu  ^ 

li'Wl.,  |y(n)U  ■  "so 

we  obtain 


^  d  ^,-EnlkT  _  (^1  _  QeC^) 

^0  "so  [>"(0)1.,  go 


(2.25) 


(2.26) 


where  use  has  been  made  of  Equation  2.9.  Rearrangement  of  Equation  2.26  leads  to 
the  result 


^  =  91^1  _  g  .  (2.27) 

So  n=.l  So 

For  iVj,  the  lowest  electronically-excited  state  (A^E^)  lies  about  6  eV  above 
the  ground  (see  Table  2.1)  and  so,  for  temperatures  ^  10,000  °K,  the  Boltzmann 
factors  in  the  summation  term  on  the  right  side  of  Equation  2.27  will  be  very  small 
(<  10~^).  Moreover,  for  all  cases  considered,  whether  for  T,  =  T  where,  it  turns  out, 
d„  <  1  or,  for  Tt  ^  T  where  d„  >  1,  the  Boltzmann  factors  are  sufficiently  small 
to  render  the  summation  term  in  Equation  2.27  small  compared  to  the  first  term. 
Therefore,  to  good  approximation. 


So 


(2.28) 


where,  tor  Nj,  go  =  1. 

Finally,  consider  the  electronic  partition  function,  Qf^*(T)  . 


Electronic  Partition  Function  for  Nj 

We  have  used,  in  definition  2.10,  the  energies  and  statistical  weights  for  the 
first  16  states  of  Wj  (up  to  and  including  6^n„)  (Reference  2.23)  to  compute  values 
for  Q^'‘(T)  for  temperatures  up  to  35,000  “K.  The  results  have  been  fitted,  to  within 
2  percent,  by  the  expression 
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gf’(r)  «  1  +  3.256  X  .  (2.29) 

For  temperatures  up  to  at  least  80(X)  “K,  Equations  2.28  and  2.29  show  that  «  1 
is  a  good  approximation. 
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SECTION  3 


DELTA  MODEL  FOR 


STATES  AND  BAND  SYSTEMS  INCLUDED 

Table  3.1  shows  the  electronic  states  of  included  in  our  model  as  well  as 
their  energies  (Ref.  2.1)  and  statistical  weights.  The  energies  are  given  with  respect 
to  both  the  ground  state  of  TVj  and  the  ground  state  of  ■ 

Table  3.1.  States  of  N^  included. 


State 

(v=0) 

Energy  Relative  to  Ni(X^Eg) 
(cm'^eV) 

Energy  Relative  to  N^(X^S^) 
(cm“^/eV) 

Statistical  Weight 

g 

12,5667.5/15.58 

2  ! 

13,4683.9/16.70 

4  ' 

15,1233.5/18.75 

2 

19,0209.5/23.58 

2  ' 

1 

Table  3.2  shows  the  three  band  systems  included.  Each  is  a  permitted  transi¬ 
tion.  Also  shown  in  Table  3.2  are  the  (0-0)  band-origin  wavelengths  and  wavenumbers 
of  the  systems  and  the  approximate  wavelength  range  over  which  each  system  radiates 
(Ref.  2.1). 


EXCITATION /DEEXCITATION  MECHANISMS 


The  mechanisms  included  for  excitation  and  quenching  of  the  electronic 
states  of  Nj  zu’e  electron  and  heavy-particle  collisions,  spontaneous  radiation  and 
sunlight  absorption.  Excitation  of  the  excited  states  of  N2  is  assumed  to  occur  only 
by  collisions  on,  and  sunlight  absorption  by,  NjiX^'L'^)  ions.  Excitation  by  collisions 
on  N2(X^T.j)  is  ignored  because,  at  the  temperatures  of  interest,  the  excitation  rate  is 
very  small  compared  to  that  by  collisions  on  N^(X^E^)  ions.  For  example,  even  for  a 
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Table  3.2.  Band  systems  of  N]  included. 


Bamd  System 

Kn 

f'oo 

(cm-i) 

Wavelength  Region 

(A) 

A^n„-X*E;  (Meinel) 

11,090.9 

9,016.4 

17,700-5,500 

B®E:-X*E;  (1st  Neg.) 

3,911.4 

25,566.0 

5,870-2,860 

C®E:-X®E:  (2nd  Neg.) 

1,549.4 

64,542.0 

2,230-1,270 

temperature  as  high  as  30,000  °K,  we  estimate  that  the  excitation  rate  of  iV2"^(A^nu) 
by  e  +  )  will  dominate  that  by  <  +  N2{X^E* )  for  fractional  ionization  as 

small  as  0.2  percent.  The  same  will  be  true,  to  am  even  greater  extent,  for  excitation 
of  S*E+  and 

Table  3.3  shows  the  specific  reactions  included  in  our  model  and  sug¬ 
gested  values  for  the  rate  coefficients,  the  Einstein  A  coefficients  and  the  pB  factors 
for  sunlight. 

Laboratory  data  on  excitation  of  the  electronic  states  of  by  electron 
impact  on  N2  (-X’^E^ )  appear  to  be  very  limited.  Most  data  refer  to  excitation  by  im¬ 
pact  on  JV2(X^E^).  However,  the  absolute  cross-section  measurements  of  Crandall  and 
Kauppila  (Ref.  3.6),  for  excitation  of  by  electron  impact  on  N^(X^Sj), 

are  within  a  factor  of  about  2  at  threshold  (3.17  eV)  and  within  20  percent  at  the 
highest  energy  (91  eV)  of  those  derived  from  the  Gaunt-factor  formula  of  Seaton  (Ref. 
3.1)  for  dipole-allowed  transitions.  We  have,  therefore,  chosen  to  use  the  Seaton  for¬ 
mula  for  the  cross  sections,  <r,  to  calculate  the  rate  coefficients  for  Reactions  3-1  to 
3-3  in  Table  3.3.  When  the  quantity  va  is  averaged  over  a  Maxwellian  velocity  (u) 
distribution,  the  rate  coefficient  for  excitation  of  state  n  from  state  m  (n  >  m)  can 
be  written  (Ref.  3.2)  as 


k„r,  =  1.57  X  (cm®  sec"*)  (3.1) 

where  the  electron  temperature,  T*,  and  excitation  energy,  =  F„  —  B„,  are  in  eV', 
G(EfxlT,)  is  an  effective  Gaunt  factor  and  /„„  is  the  absorption  oscillator  strength. 

We  find  that  the  values  for  G  tabulated  in  Ref.  3.2  for  ionic  species  can  be 
represented,  to  good  approximation  for  rt<5Ee„  by  the  formula 
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Gi{E„IT,)  =  0.2  +  0.356e"*"*^"/^‘  .  (3.2) 

Furthermore,  utilizing  the  relationship  between  oscillator  strength  and  Einstein  A 
coefficient  (Ref.  3.8) 


/mn  — 


gn  An, 


and  expressing  Tj  in  "if,  we  can  write  Equation  3.1  as 


(3.3) 


kr.n  =  (o.2  +  0.356c-=*“”«^«/^-)  (3.4) 

where  Ej,  is  in  eV.  From  detailed  balance,  the  deexcitation  rate  coefficient,  knm,  is 
given  by 


k 


nm 


AKCA.SE../T. 

^  ^  ^mn 

?n 


(3.5) 


so  that 


=  +  (cm’ sec-)  .  (3.6, 

^  e  ^ex  *^nm 

Given  the  form  of  Expression  3.5,  we  write  the  deexcitation  rate  coefficients 
3-1  through  3-3  (Table  3.3)  as 


Jk  =  a  r*(d  4-  (3.7) 

and  the  corresponding  excitation  rates  as 

k'  =  i=.e-ii604.5£../rejt  =  e-c'/T.jt  .  (3.8) 
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The  parameter  values  in  Equations  3.7  and  3.8,  for  Reactions  3-1  to  3-3,  follow  from 
Equation  3.6  and  the  energy,  wavenumber  and  Einstein  A  values  shown  in  Tables  3.1, 
3.2  and  3.3. 

It  is  difficult  to  assign  a  temperature  range  for  validity  of  the  rate  coefficients 
for  the  first  three  reactions  in  Table  3.3.  However,  to  the  extent  that  the  theoretical 
cross  sections  are  reasonably  good,  especially  for  electron  energies  <  15  eV,  the  rate 
coefficients  should  be  valid  over  a  wide  temperature  range,  including  the  one  shown 
in  Table  3.3. 

For  quenching  by  heavy  particles  (Reactions  3-4  through  3-11  in  Table  3.3), 
we  have  found  very  little  data  in  the  literature.  The  limited  data  for  the  A^Ilu  and 
states,  at  near  room  temperatures,  is  contained  in  the  references  shown  in  Table 
3.3.  The  rate  shown  for  quenching  of  A*n„  by  N  and  O  atoms  is  an  estimate  based  on 
quenching  by  He  atoms.  The  quenching  rates  shown  for  the  state  are  assumed 

to  be  the  same  as  those  for  the  state. 

The  Einstein  A  coefficients  in  Table  3.3  for  the  Meinel  and  First  Negative 
Bands  (Reactions  3-12  and  3-13)  are  taken  to  be  the  sum  of  those  connecting  the 
u'  =  0  vibrational  state  to  the  collection  of  vibrational  states  v"  of  .  For  the 
Second  Negative  Bands,  the  Einstein  values  are  less  certain,  but  the  value  adopted 
should  be  good  to  within  a  factor  of  about  4.  The  pB  factors  for  sunlight  are  based 
on  Equation  2.5  using  the  parameter  values  in  Tables  3.1  and  3.2. 


DELTA-FACTOR  EQUATIONS  FOR 

Table  3.4  shows  the  delta-factor  notation  used  for  the  states  of  N2  • 
Table  3.4.  Delta-factor  notation  for  states  of  N2 . 


State 

1^1 

m 

mm 

From  the  set  of  14  reactions  in  Table  3.3,  we  see  that  each  of  the  three 
excited  states  it,  coupled  only  to  the  ground  state  of  N^-  Determination  of  their  delta 
factors,  therefore,  is  straightforward.  It  follows  a  procedure  similar  to  that  described 
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in  the  previous  section  for  in  which  we  solve  the  steady-state  rate  equations  for 
the  excited-state  populations  in  terms  of  the  ground-state  population.  The  results  are 
as  follows: 


^  ^^I297r7/T 

k\Nt  +  k^fifi[M]  -h  An 


where,  for  the  ground-state  delta  factor,  we  use  the  approximation  (compare  Equation 
2.28) 


3s^{X) 


(3.10) 


Finally,  from  the  electronic  partition-function  definition  (2.10),  using  the 
ground  and  first  four  excited  states  £>*11,  and  C*E;^,  we  arrive  at  the 

result 


as  +  +  +  .  (3.11) 

From  this  expression,  we  cam  infer  that  begins  to  depart  significantly  from  unity 
only  above  about  5000  ’’K. 
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SECTION  4 


DELTA  MODEL  FOR  Oj 


STATES  AND  BAND  SYSTEMS  INCLUDED 


The  O2  molecule  is  a  weak  light  emitter  because  transitions  between  most 
of  its  low-lying  states  are  strongly  forbidden.  The  exception  is  the 
transition  that  gives  rise  to  the  Schumann-Runge  (S.R.)  bands  that  span  the  region 
5350  -1750  A  (Ref.  4.1). 

These  are  the  only  bands  considered  here.  Because  only  the  one  electronically- 
excited  state  of  O2  is  involved,  we  expand  our  delta-factor  determinations  to  include 
vibrational  states  u'  =  0  to  20  of  the  state. 

Table  4.1  shows  some  parameter  values  for  the  0-0  band  of  the  S.R.  system. 
Energies  for  vibrational  states  of  up  to  v'  =  20  and  wavenumbers  for  associated 
band  systems  are  given  in  terms  of  equations  cited  below. 

Table  4.1.  State  and  band-system  parameters  for  O;. 


State 

(v=0) 

Energy 

(cm“^/eV) 

Statistical 

Weight 

g 

Band  System 

•^00 

(A) 

^00 

(cm-^) 

Wavelength 

Region 

(A) 

~X^ 

3 

(S.R.) 

2026.01 

49358.15 

5350-1750 

49358.17/6.12 

3 

EXCITATION  /DEEXCITATION  MECHANISMS 


The  mechanisms  included  for  excitation  and  quenching  of  the  B^Y.~  state 
are  electron  collisions,  spontaneous  radiation  and  sunlight  absorption.  We  have  found 
no  data  in  the  literature  on  quenching  by  heavy  particles  and  so  exclude  the  process 
from  further  consideration.  However,  with  a  radiative  lifetime  of  typically  about  40 


ns,  one  would  not  expect  such  collisions  to  be  effective  in  quenching  above  about 

40  km  in  the  atmosphere. 

The  electron  excitation /deexcitation  reactions  included  in  our  delta  model 
are 


e  +  O,(X%-;t;"  =  0-23)  e  +  u' =  0  -  20)  (4.1) 

and  the  radiative  processes  are 

=  0-20)  v"  =  0  -  23)  + /it/(5.i2.)  (4.2) 

where  and  k^iv"  denote  rate  coefficients  for  excitation  and  deexcitation,  respec¬ 
tively,  and  j4„»uh  zmd  pBv»v>  the  usual  Einstein  A  coefficients  amd  pB  factors.  Note 
that  we  include  excitation/deexcitation  of  21  vibrational  states  of  and  24  vi¬ 
brational  states  of  the  ground  state  The  electron  and  radiative  processes  will 

now  be  considered  separately. 


Electron  Excitation/Deexcitation  Rate  Coefficients 

There  is  a  dearth  of  data  in  the  literature  on  the  electron-excitation  process 
4.1  involving  bound  vibrational  states  of  the  electronic  state  B^E~.  Most  cross-section 
measurements  and  calculations  (see,  for  example,  References  4.2  and  4.3,  respectively) 
include  the  absorption  continuum  which  au:counts  for  nearly  all  the  transition  prob¬ 
ability.  For  example,  the  total  electronic  absorption  oscillator  strength,  fu,  for  the 
X  -*  B  transition  is  about  0.162  (Ref.  4.4)  whereas,  for  transitions  to  bound  states, 
Hv'  fov'  ^  0.00034  (Ref.  4.5).  The  vibrational  states  of  B^E~  above  v'  =  2  are  highly 
predissociated.  The  only  data  on  electron-impact  cross  sections  to  bound  states  of 
B^E~  we  can  find  are  some  theoretical  results  of  Garrett  et  al.  (Ref.  4.6)  for  electron 
energies  up  to  25  eV.  They  calculated  cross  sections  for  the  transitions  v  — ►  v' 
for  v"  =  0, 1  and  2.  We  have  compared  these  cross  sections  with  those  derived  from 
the  Seaton  formula  (Ref.  3.1)  for  dipole-allowed  transitions  and  find  that  the  two 
results  agree  to  within  a  factor  of  about  2.  For  this  reason  we  will  use  the  Seaton 
cross  sections  and  the  excitation  rate  coefficients  for  Reactions  4.1  derived  from  them. 

The  rate  coefficients  based  on  the  Seaton  cross  sections  are  represented  by 
Equation  3.1.  In  the  present  case,  the  Gaunt  fawjtor,  G,  appeairing  in  Equation  3.1  is 
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that  for  neutral  species.  From  the  tabulated  values  for  G  given  in  Reference  3.2  we 
find,  to  good  approximation  for  T*  <  the  relation 


Gn{E„IT,)  =  (4.3) 

where  Tg  and  Eg^  aire  in  eV.  Equations  3.1  and  4.3,  with  Tg  in  “K  and  h  !/„/„»  in 
cm"^  lead  to  the  result 


k,n„.  =  0.754  477  «« (cm^  sec'^)  (4.4) 


-3/2 

t'o'u" 


where  are  the  band-origin  wavenumbers  of  the  S.R.  bands  and,  numerically,  hc/kt, 
=  1.4388  cm(‘’K). 

The  band  wavenumbers  are  determined  from  the  relation  (Ref.  2.2) 


(4.5) 


where  the  energies  (cm"^)  of  the  upper  and  lower  states  are  given,  respectively,  by 


^Ug  +  u'^iv'  +  1/2)  -  +  1/2)*  +  u;:y:(v'  +  1/2)* 


(4.6) 


and 


i>(X*E7; u")  =  u'lW'  +  1/2)  -  JXiv”  +  1/2)*  +  +  1/2)*  +  +  1/2)*  . 

(4.7) 

Values  (cm~')  that  we  adopt  for  the  electronic  and  vibrational  parameters  in  Equa¬ 
tions  4.6  and  4.7  are  (Ref.  2.23): 

Og  =  49793.28;  =  709.31;  u;'X  =  10.65;  u[y[  =  -0.139 

u'l  =  1580.193;  =  11.981;  =  0.04747;  =  -0.001273  . 
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As  for  the  absorption  oscillator  strengths  for  the  S.R.  bands  in  Equation 
4.4,  we  use  the  relation  (Ref.  4.4) 

fv“v'  —  2.36  X  lQ~^q^„nu^^„u  (4.8) 

where  q„i„»  is  the  Franck-Condon  factor  connecting  the  upper  and  lower  states.  For 
these  Framck-Condon  factors,  we  adopt  the  results  of  Harris  et  al.  (Ref.  4.5)  that 
provide  values  for  v'  up  to  20  and  u"  up  to  23. 


Einstein  A  Coefficients  and  pB  Factors 

The  Einstein  A  coefficient  connecting  vibrational  state  v'  of  with 

vibrational  state  v"  of  is  related  to  the  band  wavenumber  and  oscillator  strength 
through  Equation  3.3.  Numerically,  this  relationship  can  be  written 

=  0.667  /«<<«»  .  (4-9) 

Values  of  Au'p"  for  the  S.R.  bands  can  be  determined  from  Equations  4.5  through  4.9 
given  the  statistical  weight  factors  in  Table  4.1  that  yield  =  1. 

The  pB  factors  for  excitation  of  the  S.R.  bands  by  sunlight  absorption  are 
determined  from  Equation  2.5  using  the  A„i„i’  and  values  established  above. 


DELTA-FACTOR  EgUATIONS  FOR  Oj 


Given  the  rate  coefficients  and  optical  parameters  described  above,  we  are 
now  in  a  position  to  formulate  steady-state  equations  to  determine  delta  factors  for 
the  states  0}(fl®E~;  u'  =  0  —  20).  For  simplicity  of  notation,  we  will  hereafter  denote 
the  upper  states  as  02{B]v')  and  the  lower  states  as  02{X\v"). 

By  equating  the  production  and  destruction  rates  of  state  02{B\v'),  we 
arrive  at  the  steady-state  population 
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(4.10) 


[0j(5;  «')]„  = 


y^.  (N*  [Oi{X;  v”)],. 


00 

o"=0 


Here  the  rate  coefficient  for  electron  deexcitation  of  02{B\v'),  is  given  in  terms 
of  the  excitation  rate  i  nefficient  by  the  detailed-balance  relation 


k^,^„  =  inL  ^(E(B..')-B(xy'))/k,T.  (4  11) 

where  the  statistical-weight  ratio  is  unity  and  the  state  energies,  E,  are  in  ergs. 

Next,  since  the  Boltzmann  population  of  Oi{B;v')  can  be  written  as 


and  since  the  delta  factor  for  Oi{B\v')  is  defined  as 

«o.(«y)  -  io,{S;v')U  ' 

Equations  4.10  through  4.13  can  be  comDined  lo  yield  tne  result 


(4.12) 


(4.13) 


So,(By)  =  -  •  (4-14) 

£  +  A,.,.,) 


Here,  6o.j{x-,v")^  the  delta  factor  for  vibrational  state  v"  of  the  ground  electronic  state, 
is  defined  by  the  ratio 


<^03(X;u") 


[02{X;v")]„ 

[0,(X;v'')U 


(4.15) 
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This  ratio  can  be  reduced  by  writing  the  Oi(X;v''}  population  as  the  product  of  the 
total  02(X)  population,  [Oj(JY’)],  and  the  probability,  /*(«”),  that  vibrational  state  v" 
is  occupied.  Thus, 


^Oj(X;v") 


(o,(jf)U  p„(»') 


(4.16) 


We  next  assume  that,  in  steady  state,  the  vibrational  states  are  Boltzmann  populated 
so  that  P„(u")  =  Pe,(u'').  This  is  a  reasonable  assumption  in  view  of  the  fact  that 
O2  is  an  infrared-inactive  molecule  for  which  radiative  transitions  between  vibrational 
states  of  its  ground  electronic  state  are  forbidden.  Equation  4.16  now  reduces  to 


c  [02{X)\,t  _  ^  _  jQ  f  A  , 

Sojixy)  «  =  Soj(x}  =  %  (4.17) 

where  the  notation  for  the  ground-state  delta  factor  conforms  with  that  employed 
in  Sections  2  and  3  for  iVj  and  Wj",  respectively. 

For  we  again  invoke  the  approximation  2.28  even  though,  for  O2,  it  is 
less  valid  than  it  is  for  Wj.  This  is  because  the  first  two  electronically-excited  states  of 
oxygen  (a^A,  and  6^2^)  lie  only  about  1  and  1.6  eV,  respectively,  above  the  ground 
state  and  so  the  summation  terms  in  Equation  2.27  are  not  so  negligibly  small.  Thus, 
the  approximation  2.28  will  yield  an  upper  limit  to  the  true  value  of  d^’.  The  key 
ingredient  in  Equation  2.28  is  the  electronic  partition  function  for  Oj  that  can  be 
calculated  from  Equation  2.10  given  the  energy  and  statistical  weights  of  the  states. 
We  have  used  the  first  9  excited  states  of  O2  (up  to  c'®Aju  at  9.29  eV),  with  parameter 
values  given  in  Reference  2.23,  to  calculate  Qf^lgo  as  a  function  of  temperature  to 
12,000  “K.  The  results  can  be  fitted,  to  within  about  3  percent,  by  the  expression 


d?’  «  1  +  1.336  X  10"*  g-io309.i/r 

The  values  for  d^’  begin  to  depart  significantly  from  unity  only  above  about  5000  °K. 

The  delta  factors  for  Oj(fl;  v'  =  0  —  20)  can  now  be  determined  from  Equa¬ 
tions  4.14,  4.17  and  4.18  given  the  parameter  expressions  in  the  previous  subsection. 
The  summations  in  Equation  4.14  are  performed  up  to  v"  =  23  only. 
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SECTION  5 


DELTA  MODEL  FOR  NO 


STATES  AND  BAND  SYSTEMS  INCLUDED 


Table  5.1  shows  the  electronic  states  of  NO  included  in  our  model  together 
with  their  energies  (Ref.  2.3)  and  statistical  weights.  There  is  spin-orbit  splitting  of 
the  states  into  and  that  have  slightly  different  energies.  Similarly.  *A 

is  split  into  and  ^As/j.  However,  to  within  the  accuracy  of  the  present  model, 

we  can  ignore  these  fine-structure  splittings.  Where  distinction  is  made  in  Reference 
2.23,  we  have  used  the  energies  of  the  states. 


Table  5.1.  States  of  NO  included. 


State 

Energy  (v=0) 
(cm"V®V) 

Statistical  Weight 

g 

O.O/O.O 

4 

43,965.7/5.45 

2 

B^n 

45,913.6/5.69 

4 

c^n 

52,126/6.46 

4 

53,084.7/6.58 

2 

5'*A 

60,364.2/7.48 

4 

60,628.8/7.52 

2 

Table  5.2  shows  the  band  systems  included  together  with  the  (0-0)  band- 
origin  wavelengths  and  wavenumbers  of  the  systems  and  the  approximate  range  of 
wavelengths  over  which  each  system  radiates.  All  band  systems  are  electric-dipole 
allowed. 


EXCITATION /DEEXCITATION  MECHANISMS 


The  processes  included  here  for  excitation  and  quenching  of  the  excited 
electronic  states  are  electron  and  heavy-particle  collisions,  spontaneous  radiation  and 
sunlight  absorption.  Table  5.3  shows  the  specific  reactions  incorporated  in  our  model 
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Table  5.2.  Band  systems  of  NO  included. 


Band  System 

HU 

f'oo 

(cm"*) 

Wavelength  Region 

(A) 

(gamma) 

2274.5 

43,965.7 

1600-2700 

B^U-X^n  (beta) 

2178.0 

45,913.6 

1800-4000 

C*n-X*n  (delta) 

1918.4 

52,126 

1650-2600 

(Heath) 

12,254 

8160.3 

11,800-12,300 

(epsilon) 

1883.8 

53,084.7 

1600-2300 

£»*E+-A*E+  (Feast  1) 

10,966 

9119.0 

9200-11,500 

B*  *A-A^n  (beta  prime) 

1656.6 

60,364.2 

1500-2300 

£*E'*'-Ar*n  (gamma  prime) 

1649.4 

60,628.8 

1500-2000 

E*E+-A*E+  (Feast  2) 

6001.3 

16,663.1 

5970-6002 

f;*E+-c*n 

11,761 

8502.8 

11,700-12,200 

(Feast-Heath) 

13,255 

7544.1 

12,500-  13,300 

along  with  suggested  values  for  the  rate  coefficients,  Einstein  A  coefficients  and  pB 
factors  for  sunlight.  The  deexcitation  (forward)  rate  coefficients  in  Table  5.3  are 
expressed  in  the  form 


k  =  ae^^°^**  (cm®  sec"*)  (5.1) 

and  the  excitation  (reverse)  rates  are 

k' =  a!  k  (cm*  sec"*)  .  (5.2) 

The  quenching  data  available  for  our  purposes  leaves  much  to  be  desired. 
For  quenching  by  heavy  particles,  only  room-temperature  data  seem  to  be  available 
and  often  the  quenching  products  are  not  determined  in  the  experiments.  When 
these  products  are  unknown,  the  reverse  (excitation)  reactions  cannot  be  allowed  for 
properly.  Moreover,  for  certain  states,  especially  and,  to  a  lesser  degree, 

the  data  are  very  fragmentary. 
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Tabic  5.3.  Excitation/dcexcitation  mcchanums  and  assumed  rates  for 
states  of  NO. 
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Table  5.3.  Excitatioii/deexcitation  tncchanuiins  and  assumed  rates  for 
states  of  NO  (Coutiuued). 
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Electron- Impact  Excitation/Deexcitation 

For  electron-impact  excitation  of  the  states,  the  only  absolute  cross-  section 
data  we  have  found  in  the  literature  pertain  to  certain  bands  of  the  gamma  system 
(Ref.  5.17)  and  to  the  (4,0)  band  of  the  beta  prime  system  (Ref.  5.18).  To  see 
if  it  is  reasonable  for  us  to  use  Equations  3.1  and  4.3,  based  on  the  Seaton  cross 
sections  (Ref.  3.1),  to  calculate  the  excitation  rate  coefHcients,  we  compared  the 
results  obtained  from  the  Seaton  formula  for  the  7(1,0)  and  7(1,5)  bands  with  those 
calculated  using  Equation  2.3  emd  the  cross  sections  in  Reference  5.17.  In  both  cases, 
the  rate  coefficients  agree  to  within  a  factor  of  2  for  electron  temperatures  above  about 
3000  "K.  For  this  reason,  we  have  adopted  Equations  3.1  and  4.3  for  determination  of 
the  excitation  rate  coefficients  5.1  through  5.6  in  Table  5.3. 

When  Equations  3.1  and  4.3  are  combined,  we  can  write  the  electron-  exci¬ 
tation  rate  coefficients  in  the  form 


=  1.04  X 


(5.3) 


where  /  is  the  average  absorption  oscillator  strength  connecting  the  ground  state  and 
electronically-excited  state,  E„  is  the  average  excitation  energy  (eV)  and  Te  is  the 
electron  temperature  (®K).  We  define  the  average  oscillator  strength  and  excitation 
energy  by  the  relations 


/ 


E 


es 


oo 


y"*.  /o«'9w'o 

(5.4) 

»'=0 

i>/8065.48 

(5.5) 

OO 

£'v'o?v'o  ) 

«'=0 

(5.6) 

where  g^'o  is  the  Franck-Condon  factor  (Ref.  5.19)  connecting  the  v"  =  0  vibrational 
state  of  the  ground  electronic  state  with  the  v'  vibrational  state  of  the  excited  elec¬ 
tronic  state  and  Equation  5.5  is  the  relation  between  the  excitation  energy,  in  eV 
units  and  P  in  cm"'.  To  be  more  precise,  the  sums  in  Equations  5.4  and  5.6  should 
also  be  performed  over  v"  with  appropriate  weighting  of  the  ground  vibrational-state 
populations.  However,  measured  oscillator  strengths  reported  in  the  literature  refer 
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mainly  to  fov<-  Moreover,  the  values  of  /o„i  we  have  found  are  usually  restricted  to 
v'<5.  Therefore,  the  actual  Franck-Condon  values  used  in  Equations  5.4  and  5.6  axe 
normalized  ones  that  sum  to  unity  over  the  range  u'  =  0  to  rather  than  over 
the  range  v'  =  0  to  oo.  The  first  six  entries  in  Table  5.4  show  the  values  for  /  and 


Table  5.4.  Average  excitation  energy  and  electronic-absorption  oscillator 
strengths  for  NO  transitions. 


- 1 

Transition 

B 

WSm 

HB 

/ 

Reference 

5.612 

4.526  (4) 

6.08  (-4) 

5.20,  5.21 

B^n-x^n 

6.174 

4.980  (4) 

1.40  (-5) 

5.22,  5.23 

c*n-x*n 

6.493 

5.237  (4) 

5.50  (-3) 

5.24 

7.005 

5.650  (4) 

4.30  (-3) 

5.25 

B'  ^A-X^n 

7.923 

6.390  (4) 

1.10  (-3) 

5.18,  5.26 

E^E^-X^n 

7.985 

6.440  (4) 

1.86  (-3) 

5.15,  5.19 

- 

8.160  (3) 

4.55  (-1) 

5.27  to  5.31 

- 

9.119  (3) 

2.40  (-1) 

5.27,  5.31 

- 

1.663  (4) 

5.00  (-3) 

5.31 

- 

8.503  (3) 

3.00  (-3) 

5.31 

- 

7.544  (3) 

4.10  (-1) 

5.31 

Eex  adopted  for  the  various  transitions  and  used  in  Equation  5-3  to  compute  the  ex¬ 
citation  rate  coefficients.  The  remaining  entries  in  Table  5.4,  corresponding  to  “hot 
band”  transitions,  give  values  for  u  =  i/oo  from  Table  5.2  and,  for  /,  values  averaged 
over  those  reported  in  the  indicated  references. 


Heavy-Particle- Impact  Excitation /Deexcitation 

Re2M:tions  5-7  through  5-16  pertain  to  collisions  with  heavy  particles.  With 
the  exception  of  Reactions  5-8,  5-9  and  the  reverse  of  5-11,  they  all  refer  to  quenching 
of  the  excited  states.  As  indicated  above,  most  measurements  have  been  made  at 
or  near  room  temperature  and  the  products  of  the  reactions  are  usually  not  known. 
Therefore,  application  of  the  principle  of  detailed  balance  to  the  quenching  rate  coef¬ 
ficients  to  obtain  the  excitation  rate  coefficients,  under  the  assumption  that  the  states 
are  quenched  to  the  ground,  is  hazardous.  Consequently,  we  have  omitted  the  reverse 
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(excitation)  processes  except  in  the  case  of  Reaction  5-11.  These  omissions,  together 
with  the  adoption  of  temperature-independent  rate  coefficients  for  heavy-particle  col¬ 
lisions,  may  be  significant,  particularly  at  high  temperatures.  However,  the  effect 
of  these  approximations  is  somewhat  mitigated  by  the  fact  that  the  electron  exci¬ 
tation/deexcitation  collisions,  described  above,  tend  to  dominate  the  heavy-particle 
collisions  at  high  temperatures. 

The  measured  rate  coefficients  are  often  sensitive  to  the  particular  vibra¬ 
tional  state  considered.  In  Table  5.3,  the  assigned  values  represent  averages  over 
those  reported  in  the  references  shown  in  the  last  column. 


Radiation  Parameters 

Entries  5-17  through  5-27  in  Table  5.3  show  the  assumed  Einstein  A  coeffi¬ 
cients  and  the  associated  pB  factors  for  sunlight  absorption.  The  A  coefficients  were 
calculated  from  Equation  3.3  using  the  0  and  /  values  in  Table  5.4  and  the  pB  values 
were  calculated  from  Equation  2.5. 


DELTA-FACTOR  EQUATIONS  FOR  NO 

Table  5.5  shows  the  delta-factor  notation  used  for  the  states  of  NO. 
Table  5.5.  Delta- factor  notation  for  states  of  NO. 


State 

Delta 

State 

Delta 

KM 

~~d^ 

B'^A 

dr 

dr 

c*n 

PM 

From  Table  5.3  it  is  clear  that,  in  our  scheme,  states  A,  C,  D  and  E  are 
coupled  to  each  other  and  to  the  ground,  whereas  states  B  emd  B'  are  uncoupled 
except  to  the  ground.  In  the  following  discussion  we  treat  first  the  uncoupled  states 
and  then  the  coupled  states. 
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States  5*n  and  S'* A 


Determination  of  the  delta  factors  for  the  S*n  and  S'* A  states  follows  a 
procedure  similar  to  that  described  in  Section  2  for  We  solve  the  steady-state  rate 
equations  for  the  excited-  state  populations  in  terms  of  the  ground-state  population. 
The  procedure  is  straightforward;  only  the  results  are  presented  here.  They  are  as 
follows. 


.j,o  ^  -h pSta)e»“*/^  -h  3 

*  kzNe  +  kiolM]  -f-  Ais 

.so  ^ 

*  k^Nf  -)-  fcis[iVj]  i-  Ail 

where  is  given  by  Equation  2.22,  by  Equation  2.29  and 


kio[M\  =  kio{Ni)[Ni\  -h  A:io(0,)[0,]  kio{NO)\NO]  -h  kioiHiO)[H20\  .  (5.9) 

The  parameter  numbering  in  Equations  5.7  through  5.9  follows  that  in  Table  5.3. 

For  the  ground-state  delta  factor,  we  invoke  an  approximation  analo¬ 
gous  to  2.28  emd  write 


ss  Q^°[T)lgso[xm)  •  (5.10) 

Given  the  partition-function  definition  2.10,  the  average  state  energies  (cm“*)  in  Table 
5.4  and  statistical  weights  in  Table  5.1,  Equation  5.10  becomes 


d^'^’  «  1  +  0.5e-«®“^/^  +  e-n65:/T  ^  ^-75349/r  ^  q  5^-91291/7  (5 

where  the  summation  in  Equation  2.10  is  taken  over  the  first  5  states.  From  Equation 
5.11,  we  find  that  begins  to  depart  significantly  from  unity  only  for  T  >  15, 000‘’K. 
Thus,  for  most  applications  of  interest,  ss  1  is  an  excellent  approximation. 
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Statcc  .4*E+,  C*n,  £>*S+  and 

As  mentioned  above,  these  four  states  zu^e  coupled  to  each  other  as  well  as 
to  the  ground  and  values  for  their  delta  factors  require  simultaneous  solution  of  four 
linear  steady-state  equations  that  are  easily  derived  from  the  reaction  set  in  Table 
5.3.  In  particular,  we  find  that  the  delta  factors  ■,  and  d^°  are  related 

to  dQ°  and  to  each  other  by  the  set  of  equations 


hd^^  -b,d!^0  ^  (5  12) 

b^d^^ -byd^<^ +0.d^° +  b,id^°  =  -b^d^°  (5.13) 

biod^^  +  O.d^^  —  bnd^^  +  bi2d^^  =  —bisd^^  (5.14) 

bud^*^  +  bisd^^  +  bisd^^  —  bnd^^  =  —bi^d^^  (5.15) 

where  the  coefficients  are  defined  as  follows: 


bi  —  kiNe  +  fcrlAf]  -I-  A\j  +  +  pBh  +  pB^^ 

62  =  2(fci2(iV2j  + 

63  = 

64  = 

As  =  2[(fci'7V, +  ^5i7)c®®“®^^  +  3J(:8(iV2|di^’e-®^®^/^/Qf»] 

As  =  pB^z 

b^  =  2(/c3iV<  +  k\i\N2\dQ‘*  jQ^'*  -h  fci2[iV2]  -t-  ki3[M]  +  A19  +  A23  +  pB^f^t 
As  =  A26e-”®“/^ 

Ag  =  2\%N,  +  pfli9)e«®“«/^  +  3  ik'nIiVjK’  IQ^^\ 

bio  =  pBn 

bii  —  {kiNf  +  ki4[M]  +  A20  +  A24  +  pB27)i 

6„  =  A27e-'”®/^ 
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6,3  =  + 

614  =  pB-ii 

615  = 

4.6  = 

i>17  =  (^6-^e  +  +  M7  +  -^56  +  -Aj?) 

6i8  =  2(A:;jV, +  p5j2)e®®“^^^  . 

By  reference  to  Table  5.3,  the  quantities  ki[M\,  A:i3[M]  and  kn[M\  appearing  above 
are  defined  as: 

krlM]  =  k^{Ni)[N2\  +  kr{02)[02\  +  kj(NO)[NO]  +  kj{C02)[C02]  +  ferCiTjO)  [ffiOj 
kn[M]  =  k^2{[N2]  +  [NO]) 

ki4\M]  =  fci4(iV,)[iVj]  +  fci4(02)(02]  +  A:i4(iV0)[iV0]  +  fci4(C02)[C0j]  . 

The  delta-factor  solutions  to  Equations  5.12  through  5.15  are: 

di^  =  {^5(69(611617  —  6ij6i6)  +  6g(6i36i6  +  611618)]  +  615(63(69612  —  68613) 

+611(6469  —  6568)]  +  67(65(611617  —  612616)  +  613(64616  +  63617) 

+618(63612  +  64611)]}  fD  (5.16) 

ds^  =  {66(65(611617  —  612615)  +  613(63617  +  64615)  +  618(63612  +  64611)1  -t-  69(61(611617 

— 612615)  —  610(63617  +  64615)  —  614(63612  +  64611)1  +  68(61(613615  -I-  611613) 

+610(65615  —  63613)  +  614(65611  +  63613)]}  (5-17) 

d^^  =  {610(62(69617  +  63618)  +  67(65617  +  64613)  -i-  615(6469  —  6563)] 

—613(61(63615  —  67617)  +  65(67617  +  64615)  -t-  614(6763  +  6467)] 

+612(61(67618  +  69615)  —  65(62613  —  65615)  +  614(6769  -i-  6567)]}  (5.18) 

=  {65(615(63613  +  65611)  —  62(611613  +  613615)]  +  67(61(611618  +  613615) 

—610(63618  —  65615)  +  614(63613  +  65611)]  +  69(67(610615  +  611614) 

—615(63610  —  61611)]}  dq^  I D  (5.19) 
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whcxe 


D  =  biQ[bi{bsbis  —  f>7i>n)  —  +  b^b^)]  ~  ~ 

+bi(bibn  +  b^bn)  +  64(^^is  +  ^^14)]  —  ^uf^rC^i^ie  +  ^3^14) 
+  66(^3^»18  ~  ^^»16)1  • 


(5.20) 
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SECTION  6 


DELTA  MODEL  FOR  NO+ 


STATES  AND  BAND  SYSTEMS  INCLUDED 


The  major  band  system  of  NO'^,  and  the  only  one  considered  here,  originates 
from  the  transition  that  gives  rise  to  the  Baer-Miescher  (B.M.)  bands  in 

the  vacuum  UV.  Because,  as  for  Oj,  only  one  electronically-excited  state  of  iVO'*' 
is  involved,  we  increase  our  delta-factor  determinations  to  include  vibrational  states 
u'  =  0  to  10  of  the  state. 

Table  6.1  shows  some  parameter  values  for  the  0-0  band  of  the  B.M.  system. 
Energies  for  vibrational  states  of  ^4^11  up  to  v'  =  10,  and  wavenumbers  for  associated 
band  systems,  are  given  by  equations  shown  later  in  this  section. 

Table  6.1.  States  and  band-system  parameters  for  NO*^. 


State 

(v=0) 

Energy 

(cm-VeV) 

Statistical 

Weight 

g 

Band  System 

Aoo 

(A) 

Voo 

(cm-^) 

Vavelength 

Region 

(A) 

O.O/O.O 

1 

(B.M.) 

1368.3 

73083.5 

1200-1700 

73083.5/9.06 

2 

EXCITATION/DEEXCITATION  MECHANISMS 

The  mechanisms  included  here  for  excitation  and  deexcitation  of  the  A^IT 
state  are  electron  collisions,  spontaneous  radiation  and  sunlight  absorption.  We  have 
found  no  data  in  the  literature  on  quenching  by  heavy  particles  and,  consequently, 
ignore  the  process  here.  The  fact  that  thermal  emission  from  NO'^  will  be  signifi¬ 
cant  only  in  regions  of  high  NO'*'  concentration  and,  hence,  of  high  electron  density, 
together  with  the  fact  that  electron  excitation/deexcitation  by  electrons  is  relatively 
more  efficient  than  that  by  heavy  particles,  should  be  ample  justification  for  omission 
of  the  latter  process. 
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are 


The  electron  excitation  and  quenching  reactions  included  in  this  delta  model 


e  +  NO^{X^i:^;v-  =  0-7)  e  +  jVrC>+(^'n;  t>' =  0  -  10)  (6.1) 

and  the  radiative  processes  included  are 

A^O+(/l‘n;v'  =  0-10)  iVO+(A’‘E+;v''  =  0-7)+M5.M.)  (6.2) 

where,  again,  k„ii„i  and  denote  the  rate  coefficients  for  electron-impact  excitation 
and  deexcitation,  respectively,  and  A„>„»  and  are  the  Einstein  A  coefficients 

eind  pB  factors,  respectively.  Note  that  Reactions  6.1  and  6.2  connect  8  vibrational 
states  of  the  ground  state,  with  11  vibrational  states  of  A^II.  As  addressed 

again  below,  this  limitation  in  the  number  of  vibrational  states  v'  and  v"  considered 
corresponds  to  the  size  limitation  of  the  Franck-  Condon-factor  array  we  have  found 
in  the  literature  for  the  Baer-Miescher  system. 


Electron  ExcItation/Deexcitation  Rate  Coefficients 

In  the  absence  of  experimental  data  on  cross  sections  or  rate  coefficients 
for  Reactions  6.1,  we  again  resort  to  use  of  the  Seaton  cross  sections  embodied  in 
Equation  3.6  for  the  deexcitation  rate  coefficients.  If,  in  this  equation,  we  replace  the 
excitation  energy,  Ee,(cV),  by  its  equivalent 

=  tJ„.„»/8065.479  (6.3) 

and  recognize  that,  in  our  case, 

Afxtn  —  =  Av'v"  »  (®A) 


we  arrive  at  the  result 


(0.2  +  0.356.-* 


(6.5) 


6-2 


where  r«,  the  electron  temperature,  is  in  ®K.  From  detailed  balance,  the  excitation 
rate  coefficient  is  given  by 

(6.6) 

where  the  factor  2  is  the  statistical-weight  ratio,  ffA‘n/sx‘s+>  ^6  is  Boltzmann’s 
constajit. 

The  band  wavenumbers,  are  determined  from  the  relation  (Ref.  2.2) 

=  ^(y4^n;w')  -  (6.7) 

where  the  energies  (cm“^)  of  the  upper  and  lower  states  are  given,  respectively,  by 

j/(A^n;u')  «  Pt  +w'(u'  4- 1/2)  -a;'x'(v'  +  1/2)*  +  w'y'(v'  -I-  1/2)*  (6.8) 

and 


i5(X‘E-";  u")  «  a;:'(v"  +  1/2)  -  w;'i''(v"  +  1/2)*  +  u'ly'liv-  +  1/2)*  .  (6.9) 

Values  (cm~‘)  adopted  for  the  parameters  in  Equations  6.8  and  6.9  are  (Ref.  2.23): 

P,  =  73083.5;  =  1601.93;  =  20.207;  uj[y',  =  -0.2683 

w';  =  2376.42;  =  16.262;  u'ly”  =  -0.01562  . 

Einstein  A  Coefficients  and  pB  Factors 

The  Einstein  A  coefficient  connecting  an  upper  state  n,  v'  with  a  lower  state 
m,v"  can  be  written  in  terms  of  the  Frank-Condon  factor,  ?„'«»,  as  (Ref.  2.2) 

^  64  IP;"*!*  ,  ,  , 

An.v'-my  =  I'ny-my  - 9«V'  ,  (6.10) 

3/i  g„ 
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where  P"”*  is  the  matrix  element  of  the  electric-dipole  moment  connecting  the  upper 
and  lower  electronic  states.  Thus,  for  a  given  band  system,  Equation  6.10  can  be 
recast  in  the  form 


where  K  is  independent  of  the  vibrational  quantum  numbers.  To  determi..j  the  con¬ 
stant  K  for  the  B.M.  band  system,  we  utilize  the  result  of  Hesser  (Ref.  6.1)  in  Equation 
6.11: 


Aq„ii  =  1.8  X  10^  (sec"') 


(6.12) 


together  with  the  wavenumbers  from  Equation  6.7  and  Franck-Condon  factors  of 
Nicholls  (Ref.  6.2),  to  establish  the  result 

Axi'v"  —  6.43  X  .  (6.13) 

The  pB  factors  for  sunlight  aire  again  based  on  Equation  2.5  which,  in  the 
present  case,  becomes 


pB^iiyi  =  1.08  X  10"*/luV'(c 


.2.391x10- 


_  1)- 


(6.14) 


DELTA-FACTOR  EQUATIONS  FOR  NO+ 

Formulation  of  the  equations  to  determine  delta  factors  for  the  states  iVO" 
(>l'n;u'  =  0-10)  follows  a  procedure  similar  to  that  described  in  Section  4  for  the 
states  of  Oj.  In  particular,  the  steady-state  concentration  of  state  (A'n;v'),  in  terms 
of  that  for  state  (X'S"*";  v"),  is 


f;  +  /2B„,v)  [iVO+(X;u")]„ 

lNO-{A;v%.  =  '^ - ^ -  ,  (6.15) 

«"=0 
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while  the  equilibrium  populations  are  related  in  the  usual  way  by 


[iVO+(A;t;')U  =  {gA^nlgxy^^)[NO^[X-y)Ue 


+  .kev,i,Hlk%T 


(6.16) 


Since  the  delta  factors  are  defined  by  the  relation 


(6.17) 


Equations  6.15  through  6.17  lead  to  the  result 


^JVO+(A;«') 


\  f: 

^  o»=0 _ 


[NO^iX;v")\.s 

\NO^{X-y)U^ 


00 

“h  ^tf*w**) 

u"=0 


(6.18) 


By  analogy  with  Equation  4.16  for  the  case  of  Oj,  the  ground-state  popula¬ 
tion  ratio  in  Equation  6.18  can  be  expressed  as 


^ffO+(X;v") 


lNO^(X;v'%. 

[iVO^(X;u'’)le, 


[NO^(X)]„  P.,(v”) 
[iVO^(X)U  ‘  P.,(v") 


5.(A:)<5.(t;") 


(6.19) 


where  5(A')  is  the  delta  factor  for  the  ground  electronic  state  of  NO'*'  and  5„(u")  is 
the  delta  factor  for  vibrational  state  v"  of  the  ground  electronic  state. 

Since  NO'*'  is  infrared  active,  it  is  generally  incorrect  to  assume  a  Boltzmann- 
population  distribution  for  the  vibrational  states  u”  as  we  did  in  the  case  of  Oj.  How¬ 
ever,  for  temperatures  above  about  2000  'K  or  so,  there  is  a  simple  approximation 
that  can  be  adopted.  In  particular,  it  can  be  shown  (Ref.  6.3)  that  the  vibrational 
delta  factors  for  states  v"  >  1  can  be  represented  by  the  geometric  relation 


«S„(v'')  « 


(6.20) 
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Here,  6o  and  the  vibrational  delta  factors  for  the  v"  =  0  and  1  states,  respectively, 
of  the  ground  electronic  state,  can  be  determined  from  the  infrared  Delta  model  (Refs. 
1.1,  6.3,  6.4). 

For  adopt  the  usual  approximation  (2.28)  which,  in  the  present 

case,  yields 


.  (6.21) 

Values  for  ,  the  electronic  partition  function  for  NO'*' ,  can  be  obtained  from  its 
definition  (2.10).  In  the  present  case,  by  use  of  the  lowest  four  excited  states  (Ref. 
2.23),  we  can  recast  this  definition  in  the  form 

(T)  «  1  +  +  36"®""®/^  .  (6.22) 


After  combining  Equations  6.18  through  6.21,  we  can  write  the  final  expres¬ 
sion  for  the  NO'*'(A^n;  v'  =  0  —  10)  delta  factors  as 


V0+(/t;u')  = 


„i‘=i  \  2  /  \do 


v"-l 


(6.23) 


In  practice,  because  of  limitations  in  size  of  the  Franck-Condon  array  in  Reference 
6.2,  the  sums  in  Equation  6.23  are  taJcen  only  over  values  of  u"  up  to  7. 
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SECTION  7 


SOME  NUMERICAL  EXAMPLES 


As  mentioned  earlier,  the  Delta  models  described  in  Sections  2  through  6 
have  been  computerized  with  a  main  driver  program  DELTAVU  and  separate  subrou¬ 
tines  for  each  of  the  molecules  included.  In  this  section  we  present  three  numerical 
examples  of  delta  factors  calculated  by  code  DELTAVU  based  on  NORSE-code  out¬ 
put  of  species  concentrations  and  temperatures.  In  each  exaunple,  results  are  given 
for  both  day  and  night  conditions. 

Table  7.1  shows  delta  factors  at  1  sec  in  a  cell  at  200-km  altitude,  displaced 
140  km  horizontally  from  a  MT-class  burst  at  200-km  altitude.  The  species  concen¬ 
trations  and  temperatures  input  to  DELTAVU  are  shown  at  the  top  of  the  table.  NE 
denotes  electron  density  (cm~*)  and  T  and  TE  are,  respectively,  the  neutral  zmd 
electron  temperatures.  The  delta-factor  notation  for  iVj,  and  NO  corresponds 
to  that  shown  in  Tables  2.4,  3.4  and  5.5,  respectively.  Notice  that,  in  this  case,  the 
electron  and  neutral  temperatures  are  equal. 

Two  observations  are  apparent  from  the  results  shown  in  Table  7.1.  First, 
with  the  exception  of  dj,  the  delta  factor  for  the  delta  factors  are  all  very 

small  compared  to  unity.  This  is  because  the  radiation  is  truly  collision  limited.  At 
the  high  altitude  involved,  collisional  processes,  even  at  7292  °K,  cannot  populate 
the  states  fast  enough  to  compete  with  their  rapid  depopulation  by  radiation.  This 
results  in  state  populations  that  are  orders-of-magnitude  less  than  their  Boltzmann 
values.  The  exception  occurs  for  that  is  metastable,  with  a  lifetime  of  about 

2  sec,  for  which  its  population/depopulation  rates  are  controlled  more  by  collisional 
processes  than  by  radiative  ones. 

The  second  observation  is  simply  that  the  delta  factors  for  daylight  condi¬ 
tions  are  larger  than  those  for  nighttime.  This,  of  course,  is  due  to  the  effectiveness 
of  sunlight  absorption  in  populating  the  states.  The  effect  is  particularly  dramatic 
for  the  high-lying  Rydberg  states  cs^n„  and  of  iVj  (delta  factors  dj  and  dg, 

respectively)  that  have  large  transition  probabilities  to  the  ground  and  whose  daytime 
populations  are  enhanced  about  four  orders-of-magnitude  over  those  at  night. 

Table  7.2  shows  the  results  at  1  sec  in  a  cell  at  140-km  altitude,  directly 
above  a  MT-class  burst  at  90  km.  Due  to  the  lower  altitude  involved,  the  species 
concentrations  are  typically  one-to-two  orders  of  magnitude  larger  than  those  for  the 
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case  corresponding  to  Table  7.1  although  the  electron  density  is  considerably  lower.  A 
more  significant  difference,  however,  is  that  the  electron  temperature  is  about  six  times 
Icirger  than  the  neutral  temperature.  This  difference  manifests  itself  in  an  enormous 
increase  in  the  magnitude  of  the  delta  factors.  For  exaunple,  from  Table  7.2  we  see 
that  for  the  Cg  state  of  iV*,  the  delta  factor  dw  has  the  value  2.81  x  10^®  at  night 
and  8.52  x  10®*  in  the  day.  The  reasons  for  these  very  large  numbers  have  to  do  with 
our  definition  of  the  deltas  themselves  and  with  the  fact  that,  in  the  present  case, 
most  of  the  excitation  is  produced  by  electron  impact.  The  steady-state  population  of 
the  state,  which  forma  the  numerator  of  the  deitapfactor  ratio  (compare  Equation  2.8 
for  the  case  of  is  determined  mainly  by  electrons  at  a  temperature  of  6438 

°K,  whereas  the  equilibrium  (Boltzmann)  population,  which  forms  the  denominator 
of  the  delta-factor  ratio,  is  (by  definition)  determined  by  the  neutral  temperature  of 
1046  ‘’K.  For  the  Cg  state  of  iVj,  at  an  energy  of  14.37  eV  (see  Table  2.1),  the 
equilibrium  population  in  the  present  case  is  8.9  x  10“®®  molecules  cm“*.  As  for 
the  numerator  of  the  delta-factor  ratio,  if  there  were  no  collision  limiting,  that  is, 
if  collisions  denominated  the  excitation/deexcitation  processes,  then  the  steady-state 
population  of  iVj(cs*E„)  would  equal  its  Boltzmann  population  at  6438  ®K  which, 
in  the  present  case,  has  a  value  8.44  x  10“*  molecules  cm“®.  Thus,  the  ratio  (delta 
factor)  in  this  case  would  be  9.5  x  10®^.  However,  due  to  the  large  Einstein  A  coefficient 
for  depopulation  of  the  state  by  radiation  to  the  ground  state  (Table  2.3,  Reaction 
2.26),  collision  limiting  is  very  significant.  The  steady-state  population  is  less  than 
the  Boltzmann  population  and  the  delta  factors  for  day  and  night  conditions  are 
reduced,  accordingly,  by  factors  of  9  x  10“^  and  3  x  10“*®,  respectively,  from  the  value 
of  9.5  X  10®^  to  the  values  indicated  above  and  in  Table  7.2. 

This  example,  illustrated  by  Table  7.2,  serves  to  point  out  the  importance 
of  proper  application  of  the  delta  factors  in  converting  LTE  emission  values  to  NLTE 
emission  values,  especially  when  the  electron  and  neutral  temperatures  differ.  To 
obtain  the  proper  conversion,  it  is  necessary  that  the  LTE  emission  be  calculated 
using  the  neutral  gas  temperature,  not  the  electron  temperature. 

Table  7.3  is  an  example  showing  delta  factors  calculated  from  NORSEl-code 
output  at  1  sec  following  a  MT-ciass  burst  at  600  km.  The  air  ceil  involved  here  is 
at  140-km  altitude,  displaced  horizontally  by  50  km  from  the  burst  point.  In  this 
case,  although  the  burst  and  cell  altitudes  are  quite  different  from  those  in  Table  7.1, 
the  delta  factors  are  not  greatly  different  in  the  two  cases.  The  electron  and  neutral 
temperatures  in  Table  7.3  differ  by  less  than  0.1  percent  and  so  the  situation  with 
very  large  deltas  that  we  saw  in  Table  7.2  does  not  arise  here.  For  those  states  whose 
excitation  rates  2U'e  dominated  by  neutral  collisions,  the  associated  delta  factors  in 
Table  7.3  (di,  d:,  dg  for  N2,  and  NO)  are  leu-ger  than  they  are  in  Table  7.1  where 
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the  cell  altitude  is  higher.  Other  differences  arise  due  to  the  different  temperatures  and 
electron  densities.  However,  the  results  presented  in  Table  7.3  are  a  further  example  of 
the  effect  of  collision  limiting  in  reducing  the  LTE  emission  from  permitted  electronic 
transitions  in  the  visible  and  UV  by  typically  4  to  8  orders-<^>magnitude  at  altitudes 
above  about  100  km. 


Table  7.1.  Delta  factors  at  1  sec  in  a  cell  at  200-km  altitude,  horizontally 
displaced  140  km  from  a  MT-class  burst  at  200'km  altitude. 


r«*2.50E+09  M*l.aiE+08  NO=8.00E+fl6  C02«1.50E+0«  H2O«3.1)0E«4  NU.DOE+flS  lM.OOE+09  «*1.00Em  T=7.292£-H)3  TE*7.232E+03 
DELTA  factors  FOR  ELECTRtMC  STATES  Of  N2  (OAT) 

dO  dl  d2  d3  (X  dS  d(  d7  dS  d9  310 

l.QOE^O  8.73EHI1  8.78E-0fi  a.SlE-08  a.30E-a8  8.20E-08  S.aXH18  S.88C-08  6.34EH)8  4.32E-08  4.3SE-08 

OaTA  FAaORS  FOR  ELECTRIMIC  STATES  OF  K2f  (DAY) 
d6  dl  (S  (13 

1.34E4<0  1.4SE-44  6.71EH)C  9.3X-07 

DELTA  FAaORS  FOR  02(8  tripltt  StWA  u-  ;v<0  tt  19)  (DAY) 
v*4  «*1  v«2  v*3  v*4  (/*5  (/»7  k»9  u'lD 

2.03E-OC  2.Q3EH)(  1.93C-(IC  l.TSE-OC  1.33EHI(  1.48e-4(  l.DCE-OC  1.32E-0(  1.2X-0S  1.19EH)S  1.17E-0S 

(»»11  «*12  v»13  ((*14  ((*15  \(»U  ((*17  \(»18  0*19  (»*20 

1.13E-0S  1.09C'O(  l.a<E-(K  1.07E-0S  l.SaC-OS  1.02E-06  I.OIE-Ot  l.OOC-OC  9.89C-()7  9.82E-07 

DELTA  FACTORS  FOR  ELECTR(MC  STATES  Of  NO  (DAY) 

dO  dl  d2  d3  d4  dS  d( 

l.OOEMd  2.aEH)(  l.SSE-OS  4.07E-07  3.29E-07  3.9eCH)7  1.49E-07 

OaTA  FACTORS  FOR  NOf(A  tinsltt  PI|vi4  t«  18)  (DAY) 

u*8  u*l  ((«2  u*3  u«4  ««3  ((*<  u*7  u*6  u*9  u*10 

3.1«E-fl8  l.fl3E>87  2.34e-87  3.a0E-07  3.3SE-87  3.2SW7  3.04C-07  2.«0E-fl7  1.S2E-07  1.5a-fl7  1.24E-07 


OaTA  FAaORS  FOR  EUCTRQNIC  STATES  OF  N2  (NISHT) 

d8dld2d3d4d9dSd7d8d9dll) 
l.OOEHIP  7.28E-01  5.28E-0S  8.«1E-(I8  4.38E-89  1.08E-08  5.11E-12  2,40E-12  4,3aC-10  7.97E-12  2.37E-11 

OaTA  FAaORS  FOR  EUCTRCNIC  STATES  OF  N»  (NISKT) 
dO  dl  d2  d3 

1.34E480  1.40C-84  4.iaE-8(  4,80E-87 

OaTA  factors  for  02(8  tripltt  3i9M  «•  ;v*8  tt  19)  (NISHT) 

|(*8  »*1  v*2  v*3  i(»4  i(»S  y*8  y*9  y*lO 

3.1SE-07  3.38E-07  3.24E-87  2.aE-87  1.92E-fl7  1.87E-87  I.J4C-07  1.48E-87  1.31E-07  1.20E-87  1.17EH)7 

y»ll  y»12  y*13  y»14  y*13  y*lS  y»l7  y«18  y*19  v»20 

1.09E-87  1.00EH)7  9.S8EH)8  9.74E-88  9.28EH)8  9.0S-()8  8.8ff-88  8.71E-88  8.49CH)e  8.37E-88 

DaTA  FACTORS  FOR  EUCTROHC  STATES  OF  NO  (NISHT) 

dO  dl  02  d3  04  dS  d( 

l.OOE+80  3.42EHI7  7.S9E-07  2.79E-48  1.81E-08  1.38C-88  7.1SEH)9 

OaTA  FACTORS  FOR  N»f(A  iin^ltt  PI;v*«  tt  10)  (NISHT) 
y*8  y»l  y»2  y*3  v»4  y»9  y*<  y*7  y*9  v*10 

1.03E-08  3.4tE-88  7.S3C-88  9.79C-48  1.17E-87  1.97E-87  1.00EH)7  8.f3E-e8  3.42E-48  3.22E-98  4.18E-08 
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Table  7.2.  Delta  factors  at  1  sec  in  a  cell  at  140-km  altitude,  directly 
above  a  MT-class  burst  at  90-kin  altitude. 


«*l.49£+ll  NOM.63E+08  C02«1.09e+fl7  H20««.73€+02  S»5.24£+flJ  O4.41E+10  «»5.fia+fl7  Tsl.04«'H)3  :E=«.438£t03 

OaTA  FACTORS  FOR  OICTIMKIC  STATES  OF  N2  (DAY) 

dO  jfl  d2  d3  d4  d3  dS  d7  dS  d9  dlO 

l.OOE+00  1.5SE+23  I.S3E+23  1.43E+3S  3.11E+44  7.93E+45  1.33C+46  l.«3E+4fi  7.75E+4*  fi.37E+51  0.52E+31 

OaTA  FACTORS  FOR  ELECTRONIC  STATES  OF  N2f  (DAY) 
dt  dl  d2  d3 

l.a0E4Od  3.03EH)1  2.3aE4()7  3.3DEi2C 

OaTA  FACTORS  FOR  02(8  tripltt  SIM  U'  ;v>0  to  19)  (DAY) 
vad  g»l  l/»2  V*3  uM  ll»3  ««7  «*9  v»10 

1.44E^13  2.08EM4  1.88Eil3  1.27En6  S.£7E41(  2.74(417  1.01E418  3.22E418  9.33C418  2.43C419  S.£7E4I9 

v*!!  v*12  k*13  v*14  v*19  v«l(  v*17  u>18  v>19  v*20 

1.23E42t  2.43E420  4.30E428  8.1IOC420  1.21E4a  1.79E421  2.SlE4a  3.29(43  4.09(43  4.3E43 

OaTA  FACTORS  FOR  EL£CTR(M1C  STATES  OF  NO  (DAY) 

dt  dl  d2  (13  04  «  d( 

1.00E488  7.11E417  2.9a419  8.18(420  2.33(422  1.8e(42(  1.24(426 

OaTA  FACTORS  FOR  N04(A  sinolft  Pl!«««  t(  10)  (DAY) 

v*0  v*l  <t*2  v«3  u*4  v*3  v«6  v*7  i/aS  u*9  u*10 

3.33E429  7.01E430  9,37E43  6.80E432  4.34(433  2.02(434  9.00(434  3.4«(43S  9.ia43S  3.49(436  1.03(437 


OaTA  FACTORS  FOR  ELECTRIMIC  STATES  OF  K2  (NI6KT) 

dOdl  d2d3d4d5d6d7d8d9dl0 
1.00(400  1.40E423  1.13(423  1.43(436  6.10E448  3.79(442  8.44(438  4.72(438  1.82(442  8.10(444  2.81E44S 

OaTA  FACTORS  FOR  EUCTRONIC  STATES  OF  N24  (N18HT) 
dO  dl  d2  d3 

1.00E400  1.3E-01  1.07(406  2.10(424 

OaTA  FACTORS  FOR  02(8  SIM  «•  !v«0  to  19)  (NISIT) 

v«0  v»l  «»2  o»3  »»4  v»3  7«6  y»7  y«9  y»lO 

2.33(410  3.04(411  2.32(412  1.38(413  7.77(413  3.02(414  1.00(413  3.20(413  8.91E413  2.24(416  3.04(410 

y»ll  y»12  y»l3  y»l4  y«13  y*l6  y»17  y»l8  y*19  y*20 

1.06(417  2.06(417  3.69(417  (.40(417  9.34(417  1.33(418  1.92(418  2.49(418  3.08(418  3.39(418 

OaTA  factors  FOR  ELECTRIMIC  STATES  OF  NO  (NIGHT) 

dO  dl  d2  d3  d4  d3  d( 

1.00E400  1.13(416  8.17E418  8.91E417  1.87E419  9.89(422  6.(9(422 

OaTA  FACTORS  FOR  N84(A  Pl;v4  to  10)  (NIGHT) 

y»0  y«l  y«2  y»3  y*4  y«9  y*6  y*7  y«8  y»10 

1.44(427  2.90(428  3.74(429  2.62(430  1.62(431  7.29(431  3.17(432  1.19(433  3.06(433  1.14(434  3.30E434 


7-5 


Table  7.3.  Delta  factors  at  1  sec  in  a  cell  at  140-kni  altitude,  horizontally 
displaced  50  km  from  a  MT-class  burst  at  600-km  altitude. 


M2s5. 436+10  02s3.0X+fl9  N0»3.72E+fl9  COM. 416+04  H20*5.ai6+01  M«2.776+10  O3.476+10  N6»3.806+08  1*8.9276+03  16*8.3986+03 
06LTA  FAaORS  FOR  6L6CTR0NIC  STAT6S  OF  N2  (DAY) 

dO  dl  d2  d3  d4  dS  di  d7  d8  d9  310 

1.006+00  1.646+00  2.626-05  1.086-07  4.416-09  7.806-09  1.626-09  1.606-09  1.596-09  6.8X-10  7.446-10 

OaTA  FACTORS  FOR  6L£CrR(MlC  STATFS  OF  N2+  (OAY) 
dO  dl  d2  (13 

1.486+00  1.176-03  9.366-06  4.876-06 

06LTA  FAaORS  FOR  02(8  tripitt  SIM  u-  ;v>0  t*  19)  (OAY) 
g*0  u*l  0*2  0*3  0*4  0*9  0*6  o*7  o*8  o*9  0*10 

7.436-07  7.446H)7  6.946-07  5.856-07  4.766-07  4.556-07  3.956-07  3.796-07  3.446-07  3.196-07  3.096-07 

0*11  0*12  0*13  0*14  0*15  0*16  0*17  o*18  o*19  o»20 

2.926-07  2.746-07  2.636-07  2.666-07  2.536-07  2.486-07  2.436-07  2.406-07  2.356-07  2.326-07 

OaTA  FAaORS  F(»  a6aR(MIC  STAT6S  OF  NO  (MY) 

dO  (H  d2  (13  d4  d5  d6 

1.006+00  1.176-05  3,506-05  2.406-07  4.806-08  4.556-08  1.656-08 

OarA  FACTORS  FOR  N(H(A  sin^ltt  Pl!u*0  t«  10)  (MY) 

K*0  M*1  (»»2  i<*3  v*4  g»5  g*6  g»7  g*8  g*9  g*10 

5.406-09  1.776-08  3.906-08  4,936-08  5.816-08  5.286-08  4.896-08  4.166-08  2.596-08  2.476-08  1.966-08 


OaTA  FAaORS  FOR  aSCTRtWIC  STAT6S  OF  N2  (NI6HT) 

dO  dl  d2  (13  d4  dS  d6  d7  d8  d9  dlO 

1.006+00  1.486+00  2.1X-05  1.086-07  2.296-09  6.126-09  1.676-11  8.266-12  2.i26-l0  2.316-11  8.296-11 

OaTA  FACTORS  FOR  aSCTROIIC  STAT6S  OF  N2+  (NISHT) 
dO  (U  d2  (13 

1.486+00  1.176-03  3.256-06  4.826-06 

OaTA  FACTORS  FOR  02(8  tripltt  SIM  u-  ;v4  to  19)  (NISMT) 
u*0  v*l  g«2  u*3  g*4  »*3  v*6  v*7  g*9  u*10 

1.346-07  1.446-07  1.386-07  1.116-07  8.246-08  8.006-08  6.596-08  6.376-08  5.646-08  5.156-08  5.016-08 

g*ll  v»12  v»13  g*14  w»15  v»16  g»l7  g»18  g*19  g»20 

4.696-08  4.336-08  4,136-08  4.206-08  3.996-08  3.916-08  3.826-08  3.786-08  3.676-08  3.616-08 

OaTA  FACTORS  FOR  aSCTRIMIC  STATES  OF  NO  (NISKT) 

dO  dl  d2  (13  04  dS  (16 

1.006+00  1.016-05  3.146-05  1.686-07  7.596-09  7.876-09  2.996-09 

OaTA  factors  for  N0+(A  fin^lit  PI  ;v<0  tp  10)  (N16KT) 
g*0  g»l  g»2  »»3  g»4  v»5  v*6  g*7  g»8  u*9  g*10 

3.356-09  1.286-08  2.886-08  3.706-08  4.426-08  4.066-08  3.826-08  3.286-08  2.066-08  1.996-08  1.596-08 
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SECTION  8 


COMMENTS  AND  RECOMMENDATIONS 


By  way  of  comment,  we  would  like  to  re-emphasize  certain  points  inferred 
from  the  foregoing  sections.  First,  as  in  most  models,  assumptions  and  approxima¬ 
tions  have  been  made  that  limit  the  accuracy  of  the  results  in  one  way  or  another. 
Specifically,  with  respect  to  the  wavelength  dependence  of  the  NLTE  emission  from 
a  given  band  system  of  jVj,  and  NO,  the  largest  error  probably  stems  from  our 
assumption  that  all  bands  of  the  system  have  the  same  delta  factor.  That  is,  we  have 
treated  these  «hree  molecules  as  if  they  were  vibrationless.  To  do  otherwise,  would 
involve  a  substantial  expansion  in  the  scope  of  the  present  model  because  of  the  large 
number  of  vibrational  states  of  each  electronic  state  for  which  delta  factors  would  need 
to  be  provided.  It  seems  to  us  that  an  expansion  in  scope  such  as  this  is  not  presently 
justified.  This  is  mainly  because  the  NORSE  code  is  not  designed  to  take  advantage 
of  this  level  of  detail  in  converting  LTE  emission  to  NLTE  emission.  NORSE  cannot 
distinguish  contributions  from  each  band  of  a  given  band  system  to  the  LTE  emission 
in  a  specified  wavenumber  interval.  Moreover,  the  required  excitation/deexcitation 
rate  coefficients  connecting  the  upper  and  lower  states  of  most  bands  of  each  system 
are  not  yet  known. 

How  much  error  is  introduced  by  this  single  delta-factor  assumption  for 
electronic  states  of  iVj,  and  NO  ?  It  is  hard  to  generalize,  although  some  idea  can 
perhaps  be  obtained  by  referring  to  the  results  in  Tables  7.1  through  7.3  for  the  cases 
of  Oj  and  NO'^.  For  Oj,  our  model  provides  delta  factors  for  21  vibrational  states  of 
and,  for  NO'^ ,  11  vibrational  states  of  A^II.  In  Tables  7.1  and  7.3,  where  the 
electron  and  neutral  temperatures  are  nearly  equal,  the  vziriation  in  the  deltas  over  the 
vibrational  states  of  Oj  is  at  most  a  factor  of  4  while,  for  NO'^ ,  it  is  at  most  a  factor 
of  10.  For  the  case  represented  by  Table  7.2,  however,  where  the  electron  temperature 
is  much  greater  than  the  neutral  temperature,  the  variation  in  the  deltas  for  Oj  and 
for  NO'^  is  several  orders  of  magnitude.  Therefore,  depending  on  circumstances,  it 
appears  that  the  error  resulting  from  the  assumption  that  the  deltas  for  all  bands  of 
a  given  band  system  are  the  same  can  range  from  small  to  very  large.  Large  errors 
are  likely  to  occur  only  when  the  electron  and  neutral  temperatures  differ  significantly 
from  one  2Lnother. 

Other  sources  of  error  include  assumptions  made  on  quenching  rate  coeffi¬ 
cients  for  which  no  laboratory  data  are  available.  For  electron-impact  excitation  of 
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the  states,  we  have  made  liberal  use  of  the  Seatcn  formulas  —  an  approximation  that 
is  accurate  to  within  a  factor  of  about  2  when  compared  with  the  limited  laboratory 
data  that  do  exist.  As  for  quenching  by  heavy  particles,  reliable  rate-coefficient  data 
are  harder  to  come  by.  Where  no  data  were  found,  we  omitted  the  process  entirely 
rather  than  including  it  with  a  guess  on  its  rate  coefficient.  Assessment  of  the  impact 
of  such  omissions  on  the  results  would  require  a  sensitivity  study  that  we  have  not 
had  time  to  perform. 

A  second  comment  that  needs  to  be  emphasized  here  relates  to  the  delta- 
factor  definition  used  in  our  models,  and  to  the  proper  application  of  the  deltas  in 
converting  the  computed  LTE  emission  to  NLTE  emission.  The  delta  factor  for  each 
radiating  state  considered  is  defined  as  the  ratio  between  the  actual  steady-state  pop¬ 
ulation  of  the  state  to  its  Boltzmann  population  —  where  the  latter  is  evaluated  at  the 
neutral  gas  temperature.  In  view  of  this  definition,  it  is  very  important  that  the  LTE 
emission  that  we  wish  to  convert  to  NLTE  emission  be  computed  using  the  neutral 
gas  temperature.  For  cases  where  the  neutral  and  electron  temperatures  are  nearly 
equal,  it  clearly  doesn’t  matter  which  is  used.  However,  when  the  electron  temper¬ 
ature  greatly  exceeds  the  neutral  temperature,  as  in  Table  7.2,  then  the  calculated 
NLTE  emission  may  be  orders-of-magnitude  too  large  if  the  electron  temperature, 
rather  than  the  neutral  temperature,  is  used  for  determination  of  the  LTE  emission. 

Finally,  by  way  of  recommendations  for  improvement  and  expansion  of  the 
present  model,  we  suggest  the  following  items. 

(1)  Addition  of  other  species  whose  visible/UV  radiation  in  a  nuclear 
environment  may  be  significant.  Such  species  should  include  the  atoms  O  and  N  and 
their  corresponding  ions  O'^  and  N'^.  Emission  from  the  metastable  states  of  these 
species  has  recently  been  added  to  NORSE,  but  thermal  emission  from  their  permitted 
states  may  also  be  important  in  high-  temperature  regions.  Other  species  that  should 
be  considered  for  addition  are  CO,  COj,  NOj  and  HjO.  As  for  O3,  because  of  its 
low  binding  energy  (1.05  eV),  and  hence  propensity  to  dissociate  at  even  moderate 
temperatures,  it  is  unlikely  to  be  an  important  thermal  emitter  in  the  visible/UV 
region  except,  possibly,  under  conditions  similar  to  those  corresponding  to  Table  7.2 
where  the  neutral  temperature  is  low  and  the  electron  temperature  is  high. 

(2)  Upgrading  the  rate  coefficients  is  an  obvious  area  where  improvements 
can  and  should  be  done  periodically.  In  this  context,  three  items  may  be  worth 
mentioning.  First,  the  pB  factors  for  sunlight  absorption  in  the  present  models  are 
input  as  constants,  independent  of  altitude.  It  is  not  difficult  to  allow  for  attenuation 
of  sunlight  by  the  intervening  ambient  atmosphere  by  introducing  altitude-dependent 
values  of  pB  for  each  of  the  bands  considered.  Second,  as  new  data  on  rate  coefficients 
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for  quenching  reactions  included  in  the  present  model  become  available,  it  is  a  simple 
matter  to  make  the  data  changes  in  the  appropriate  subroutines,  and  this  should  be 
done.  Beyond  that,  when  data  become  available  on  processes  not  included  in  the 
current  models,  such  as  certain  quenching  reactions  involving  heavy  particles,  the 
appropriate  models  should  be  modified  by  adding  the  processes  to  the  production- 
and  destruction-rate  terms  for  the  states  in  question  following  the  procedure  outlined 
in  the  foregoing  sections.  Third,  where  laboratory  data  are  lacking,  it  is  probably 
worthwhile  to  perform  sensitivity  studies,  using  the  codes  developed  here,  to  assess 
the  importance  of  each  reaction  on  the  resulting  delta  factors.  Such  studies  can  be 
useful  in  helping  to  order  the  priorities  of  laboratory  research  programs. 

(3)  If  future  developments  permit  the  NORSE  code  to  dbtinguish  the 
contributions  to  the  LTE  emission,  in  a  given  wavenumber  interval,  from  the  sep¬ 
arate  bands  of  each  band  system,  then  the  way  is  open  for  expanding  the  present 
model  to  include  delta  factors  for  a  whole  array  of  vibrational  states  of  each  elec¬ 
tronic  state.  This  would  provide  for  a  much  improved  wavelength  dependence  of  the 
NLTE  visible/UV  emission.  As  mentioned  eeirlier,  however,  such  a  model  would  be 
rather  unwieldy  and,  moreover,  we  would  still  be  faced  with  the  need  to  provide  exci¬ 
tation/deexcitation  rate  coefficients  for  a  large  number  of  states  for  which  laboratory 
data  might  largely  be  lacking.  Nevertheless,  for  altitudes  above  about  100  km,  where 
electron  processes  are  likely  to  dominate,  this  need  could  perhaps  be  satisfied  suffi¬ 
ciently  by  use  of  the  Seaton  approximation  for  rate  coefficients  that  we  employed  to 
advantage  in  the  present  model.  Thus,  expansion  of  the  model  to  include  many  more 
vibrational  states  is  probably  feasible  if  the  need  appears  to  warrant  it. 
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APPENDIX 


FORTRAN  LISTING  OF  PROGRAM  DELTAVU  AND 
ASSOCIATED  SUBROUTINES 


DELWU 

C  THIS  PROGiWI  COMPUTES  DELTA  FACTORS  FOR  SELECTED  ELECTRONIC  STATES  OF 
C  AIR  MOLECULES  AND  IONS.  ONLY  ELECTRIC-  DIPOLE  ALLOWED  TRANSITIONS  ARE 
C  INCLUDED.  THE  STATES  (t€  TRANSITIONS  INaUOED  CORRESPOND.  GENERALLY, 

C  TO  THOSE  IN  THE  AfWL  CODE.  THE  MOLECULAR  SPECIES  mUDED  ARE:  N2,N2f, 

C  02, NO  ANO  NOf.  WITH  THE  EXCEPTION  OF  02  NOt,  A  SINGLE  DELTA  FACTOR  IS 
C  COMPUTED  FOR  EACH  ELECTRONIC  STATE  C0NSI0£RED-Y«O  SO  THE  SW1E  DELTA 
C  factor  IS  ASSUMED  TO  APPLY  TO  ALL  VIBRATIONAL  STATES  Of  A  GIVEN  ELECTR- 
C  QNIC  STATE.  THUS,  IN  C0WERTIN6  LTE  VISia.£/UV  EMISSION  TO  NLTE  EMISSION, 
C  THE  SANE  DELTA  FACTOR  IS  APPLIED  TO  ALL  BANOS  Of  A  GIVEN  BMID  SYSTEM. 

0  THE  EXCEPTIONS  ARE  02  ANO  N(H.  FOR  02,  SEPARATE  DELTAS  ARE  COMPUTED 
C  FOR  EACH  OF  THE  VIBRATIONAL  STATES  0,1, 2,™, 20  OF  ELECTRONIC  STATE 
C  02(B  iripletSIGMA  u-)  THAT  GIVES  RISE  TO  THE  SCHUMWW-RUNGE  BANO  SYSTEM. 

C  FOR  NOf,  SEPARATE  DELTAS  ARE  COMPUTED  FOR  EACH  Of  THE  VIBRATIOW.  STATES 
C  0,1,2,  — ,10  OF  ELECTRONIC  STATE  NOf(A  singlet  PI)  THAT  GIVES  RISE  TO 
C  THE  BAER-MIESCHER  SYSTEM. 

DIMENSION  DN2(11) ,DN2P(4) ,002(21) ,DNO( 7) ,CNOP(ll) 

C  SPECIES  CONCENTRATIONS  FOR  N2,O2,NO,CO2,H2O,N,O,ELECTR0NS. 

DATA  XN2,02,XNO,C02,H20,XN,0,XNE/1.49E11,2.20E10,4.63ES,1.09E7, 

1  6.79E2,5.24E9,4.41E10,5.66E7/ 

C  NEUTRAL  ANO  ELECTRON  TEMPERATURES 

C  CAUTION  1 ! ! !  I  FOR  MACHINES  THAT  ACCEPT  NUMKRS  ONLY  SELLER  T)W4  ABOUT 
C  lfl«38,  OVERaOWS  IN  SUBROUTINE  DELN2  MAY  OCCUR  FOR  TEMPERATURES  <  1600 
C  (DEG.  K)  (KESS  DOUBLE  PRECISION  IS  USED. 

DATA  T,TE/1046.,6438./ 

0  DO  ANO  01  ARE  DELTA  FACTORS  FOR  VIBRATIONAL  STATES  0  (W)  1  OF  THE  GROUND 
C  ELECTRONIC  STATE  OF  NOf  AS  DETERMINED  BY  THE  INFRARED  DELTA  MOOa. 

DATA  00,01/2.7,4.36E-3/ 

C  DAY(K=1)/NIGHT(K=0)  INDEX. 

K=fl 

IF(K.EQ.1)WRITE(1,20)  XN2,02,XNO,C02,H20,XN,0,»C,T,TE 
IF(K.E0.1)HRITE(1,25) 

IF(K.E0.0)WRITE(1,26) 

HRITE(1,30) 

CALL  DELN2(XN2,O2,XNO,XN,O,XNE,T,TE,K,0N2) 

WRITE(1,40)  (DN2(I), 1=1,11) 

CALL  OCLN2P()iN2,O2,XN,O,XNE,T,TE,X,0N2P) 

IF(l(.EQ,nHRITE(l,45) 

IF(K.EQ.0)HR1TE(1,46) 

WRITE(1,50) 

HRITE(1,55)  (DN2P(I),I=1,4) 

CALL  Oa02(XNE,T,TE,K,D02) 

lF(K.EQ.l)miTE(l,60) 

rF(K.EQ.0)WRITE(l,61) 

WRITE(1,65) 

HRITE(1,70)  (D02(I),I=1,11) 


iJRITEa.’Z)  (002(1),  1=12, 21) 

:F(K.E0.1)WRITE(1,80) 

IF()(.£0.i))WRITE(l,81) 

WRITE(1,85) 

CALL  D£LNO(XN2,02,XNO,C02,H20,XNE,T,TE,DN2,R,»«) 

WRITE(1,90)  iDN0(l),l=l,7) 

IF(K.£Q.1)WR1TE(1,95) 

IF(K.EQ.0)WR1TE(1,96) 

WRlTEd.lOO) 

CALL  OELNOP(XNE,T,TE,K,DO,01,DNOP) 

WR1TE(1,70)(DNOP(1),1=1,11) 

20  F0RMAT(5X,'N2='1PE8.2,'  02=dPE8.2,'  N0='1PE8.2,'  C02='1PE8.2,'  H2 
10=dPE8.2,'  N='1PE8.2,'  0='1PE8.2,'  NE*'1P£8.2,'  T='1PE8.3,'  T£='l 
2PE9.3/) 

25  FORHATdlX, 'DELTA  FACTORS  FOR  ELECTRONIC  STATES  OF  N2  (DAY)') 

26  FORHATdlX, 'DELTA  FACTORS  FOR  ELECTRONIC  STATES  OF  N2  (NIGHT)') 

30  F0RHAT(8X,'  dO  dl  d2  d3  d4 

1  dS  d6  d7  d8  d9  dlO') 

40  FORMAT(9X,lldPE10.2)/) 

45  FORHATdlX, 'DELTA  FACTORS  FOR  ELECTRONIC  STATES  OF  N2+  (DAY)') 

46  FORHATdlX, 'DELTA  FACTORS  FOR  ELECTRONIC  STATES  OF  N2+  (NIGHT)') 

50  F0RHAT(8X,'  dO  dl  d2  d3') 

55  F0RHAT(9X,4dPElfl.2)/) 

60  FORHATdlX, 'DELTA  FACTORS  FOR  02(8  tripltt  SlfiHA  u-  iw»0  to  19)  (D 
lAY)') 

61  FORHATdlX, 'DELTA  FACTORS  FOR  02(8  tripUt  SIGHA  u-  |v=0  to  19)  (N 
IIGHT)') 

65  F(»1AT(8X,'  v=0  0*1  0*2  o*3  o*4 

lo=5  0=6  0*7  0=8  0*9  o*10') 

70  FORHAT(9X,lldPE10.2)/) 

72  FORHAT(9X,10(1PE10.2)/) 

75  F0RHAT(8X,'  o=ll  v*12  v=13  v=14  v*15 

lo=16  0=17  0*18  0*19  0*20') 

80  FORHATdlX, 'DELTA  FACTORS  FOR  aECTRCNIC  STATES  OF  NO  (DAY)') 

81  FORHATdlX, 'DELTA  FACTORS  FOR  aECTRONIC  STATES  OF  NO  (NIGHT)') 

35  F0RHAT(8X,'  dO  dl  d2  d3  d4 

1  d5  d6') 

90  FORHAT(3X,7(1PE10.2)/) 

95  FORHATdlX, 'DELTA  FACTORS  FOR  MOKA  singlet  PI;o»0  to  10)  (DAY)') 

96  FORHATdlX, 'OaTA  FACTORS  FOR  NO+(A  singlet  Plio*0  to  10)  (NIGHT)' 

1) 

100  F0RHAT(8X,'  o=0  o*l  o=2  o*3  o*4 

lo=5  0=6  0*7  0*8  0=9  o*10') 

END 
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iUBROUTINE  0ELN2(AN2.02,,<}<0,XN,0,M,T.n,K,0N2) 

C  THIS  SUBROUTINE  COMPUTES  DELTA  FACTORS  FOR  11  ELECTRONIC  STATES  OF  N2. 
DIMENSION  !1(29j, RP(23),0N2(11),A(19), 8(19), C(W, 0(19), F(19), 

1  A»^a0),RB(10) 

C  ELECTRONIC  PARTITION  FUNCTION  FOR  N2. 

3£N2(X)=1. +3. *EXP(-71587.2A)+6.*(EXP( -85330. 6A)+EXP( -85435. 9A)) 

1  f3.»EXP( -94748. 4A)+EXP( -97463. 3A)+2.*(EXP( -99207, OA)+D(P( -1031 

2  59. 8A))+5.*EXP( -112082. 5A)+6.*(EXP( -125319. 5A)+£XP(-128021,3A 

3  ))+10, *EXP( -134096. 2/X)+3.*£XP(-1378fl0.4A)+6.*£XP( -140432. 0A)+ 

4  EXP( -142211. 4A)+2.*EXP(-145056.2A) 

C  TVI0-800Y  RATE  COEFFICIENTS  FOR  CaLlSIOWL  OEEXCITATICN  (Ith  REACTION) 
C  ARE  IN  THE  FORM:  R(I)=A(I)*T**e(I)*EXP(-C(I)A).  THE  RATE  COEFFICIENTS 
C  FOR  THE  REVERSE  (EXCITATION)  REAQIGNS  ARE  IN  THE  FORM:  RP(I)=0(I)* 

C  EXP(-F(I)/T)*R(I).  HERE,  T  MAY  BE  EITHER  TI€  NEUTRAL  OR  THE  ELECTRON 
C  TEMPERATURE. 

'data  A72.18E-13,6.12£-11,1.59E3,5.85E4,1.78E-15,5.5E2,2.6E4,6.68E- 

1  9,1.2E-10,1.5E-13,1.£-11,2.8E-10,2.7E-11,1.1E-10,2.4E-10,4.6E-12, 

2  6.9E-ll,2.1£-ll,4.6E-7/ 

DATA  B/, 861,. 427, -2, 469, -2. 681, 1,068, -2. 469, -2. 681, -.102,. 2245, 

1  .8539, 8*0. ,-1,5/ 

DATA  C/385. 1,6759. 8, 97547. 3, 101251., -4454. 7, 92502. 4, 101251. 6, 

1  3599.9,2297.1,-2225.5,9*0./ 

DATA  0/2., 1., 2., 1., 2., 2,, 1., 6., 6., 3., 5*6., 4*3./ 

DATA  F/145054.9,14916fl.2, 149832.7,150099. ,152094.7,166375.4, 

1  166720.9,128020.1,85329.9,71586.5,2*128020.1,3*85329.9,4*71586.5/ 

C  (V)(I)  ARE  THE  EINSTEIN  'A'  COEFFICIENTS  (SEC-l)  RB(I)  ARE  THE 
C  C0RRESP0NO1N6  RH0*8  FAaORS  (SEC-1)  FOR  SIH.16HT  CONDITIONS. 

DATA  ^2*1.E8,1,4E8,1.E9,1.E8,6.2E7,4.5E7,2.E7,1.4E5,0.5/ 

DATA  R8/3.7E-8,9.3E-9,2.3E-8,7.9E-8,l.lE-8,6.6E-10,2.3E-10,9.£-2, 

1  .17,5.5E-11/ 

00  10  I»l,10 

R( I ) =A(I )*TT**e( 1 )*EXP( -C( I )/TJ) 

RP(I)=0(I)*£XP(-F(I)/TT)*R(I) 

10  CONTINUE 
DO  15  1=11,19 

R(I)=A(1)*T**B(1)*EXP(-C(I)/T) 

RP(I)=0(I)*£XP(-F(I)/T)*R(1) 

15  CONTINUE 
DO  20  1=20,29 
J=I-19 
R(I)=<VR(J) 

RP(I)=aQAT(K)*R8(J) 

20  CONTINUE 

C  DN2(1)  IS  THE  DELTA  FAaOR  FOR  TIC  GROUND  ELEaRCNIC  STATE  (X  singlet 
C  SIGMA  gf j  OF  N2,  IT  IS  APPROXIMATELY  EQUAL  TO  THE  ELECTRONIC  PARTITION 
C  FUNCTION. 

»I2(1)=0£N2(T) 

Al=(  £XP(  71586,5*(1  ,/T-l  ./H)  )*R(  10)*XNE+R(  16)*02+R(  17)*XN0+R(  18)*0 
1  +R(19)*XN)*DN2(1) 
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1 10  )tXNE+R(  16)*02+R(  1 7^*)iN0+R(  18)*(>tR(  19)*XNfR(  29) 
A3=2.*R(28)*EXP(-13743,4/T) 

Bl=iR(9)*xNE*E<P(35329.9*(l./T-l./TT)i+R(l3)*XN2+R(14)*02+R(15) 

1  ♦)(N0)*W2(1) 

Ei2=RP(28)*£XP(13743.4/'T)/2. 

B3=R(9)*XNE+R(13)*)(N2+R(14)*02+R(15)*)(N(M(28) 

e4=R(27)*EXP(-42690.2/T) 

C1=(R(8)*XNE*EXP(128020.1*(1./T-1./TT))+R(11)*)(N2+R(12)*02) 

1  *CN2(1) 

C3=RP(27)*£XP(  42690, 2/'T) 

C4=Rl 8)*XNE+R( 11 )*XN2+R( 12)*02tR( 27) 

G=A2*( B3*C4-R( 27)*RP( 27) ) -A3*82*C4 

C  THE  DELTA  FACTORS  0N2(I)  (1=2  TO  11)  REFER,  RESPECTIVELY,  TO  DELTAS  FOR 
C  THE  5  VALENCE  STATES:  A  triplet  SIGMA  u+,  B  triplet  PI  g,  C  triplet  PI  u, 
C  b  singlet  PI  j,  b'  singlet  SIGMA  u+,  FNO  THE  5  RYDBERG  STATES:  c3  singlet 
C  PI  u,  c4'  singlet  SIGMA  u+,  o3  singlet  PI  u,  c4  singlet  PI  u,  c5'  singlet 
C  SIGMA  u-f. 

CN2(  2)  =(  Al*(  B3*C4-fi(  27)*RP(  27)  )-W3*(  Bl*C4+S4itCl)  )/G 

DN2( 3) =( Al*B2*C4tA2*( B1*C4+B4*C1 ) )/G 

DN2(  4) =(  A1*B2*C3M2*(  B3*C1+B1*C3)  -R(  28)*RP(  28)*C1  )/G 

E=0, 

C  CAUTION!  III!  FOR  MACHIhCS  THAT  ACCEPT  NUMBERS  ONLY  SMALLER  THAN  ABOUT 
C  10^*38,  THE  FOLLOUING  LINES  MAY  PRODUCE  OVERFLOU  FOR  TEMPERATURES  <  1600 
C  (DEG.  K)  UNLESS  DOUBLE  PRECISION  IS  USED. 
IF(K.EQ.1)E=EXP(145054.9/T+AL0G(.5*RB(1))) 

DN2(  5) »( R(1  )*)(NE*£XP(  1 45054 , 9*(  1  ./r-1  ./TT)  )+£)/(  R(  1  )*XNE+R(  20 ) ) 

1  *0N2(1) 

E=0. 

I F( K . EQ . 1 ) E=EXP( I 491 60 . 2/T4AL0G( RB( 2) ) ) 

CN2( 6) =( R( 2)*XNE*EXP( 149160 . 2*(1 ./T-1 ./TT) )+£)/( R( 2)*XNE+R( 21 ) ) 

1  *0N2(1) 

E=0. 

I F ( R . £Q . 1 ) E=£XP( 1 49832 . 7/T+AL0G( . 5*R8( 3) ) ) 

DN2(  7) =(  R(  3)*XNE*EXP(  149832. 7*(  1  ./T-1  ./TT)  )+£)/(  R(  3)*)(NE+R(  22) ) 

1  *0N2(1) 

E=0, 

I F( K . EO. 1 ) E=EXP( 1 5O099./T+ALOG(RB( 4) ) ) 

0N2(  8)=(  R(  4)*XNE*EXP(150099. 0*(  1  ./T-l  ./TT)  )+E)/(  R(  4)*)(NE+R(  23) ) 

1  *0N2(1) 

E=0. 

: F( K . £Q . 1 ) E=£XP( 152094. 7/T+AL0G( . 5*R6( 5 ) ) ) 

DN2(  9)  =(  R(  5)*XNE*EXP(  152094, 7*(  1  ./T-1  ,/TT)  )+£)/(  R(  5)*)(NE+R(  24) ) 

1  *CN2(1) 

E=0. 

I F ( K . EQ . 1 ) £=EXP( 1 66375 . 4/T4AL0G( . 5*R6( 6) ) ) 

DN2(  10 ) =(  R(  6)*XNE*£XP(  166375. 4*(  1  ./T-1  ./H)  )+E)/(R(  6)  A)(NE+R(25) ) 

1  *0N2(1) 

E=0. 

1 F ( K . £0 . 1 ) E=£XP( 1 66720 . 9/T+AL0G( RB( 7) ) ) 
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>J2(ll)=(R(7)tXNE*E<Pabe720.9*(l./T-l.,TT))+£)/fR(7)*XNE+R(26)) 
1  ♦DN2(1) 

RETURN 

END 


SUBROI/TINE  DELN2P(XN2,02,»^,0,M,T,n,K,CN2P) 

C  THIS  SUBROUTINE  COMPUTES  DELTA  FAQORS  FOR  4  ELECTRONIC  STATES  Of  N2+. 
DIMENSION  R(14),RP(14),DN2P(4),A(11), 8(11), C(11),D(11), Gill), 

1  01(11), Fill)  ,«>(3),RB(3) 

C  ELECTRONIC  PARTITION  FUNCTION  FOR  N2+. 

QEN2P(X)=2.*(1.+2.*EXP( -12972. S8A)+EXP( -36784, 02/X)+2.*EXP(-75057 
1  .3/X)+EXP(-92862.18A)) 

C  TWO-BODY  RATE  COEFFICIENTS  FOR  COLLISIONAL  DEEXCITATION  (Ith  REACTION) 

C  ARE  IN  THE  F(»M: 

C  R( I )*A( I )*T**8( I )*(0( I )+6(I )*EXP(-C( I )/T) ) . 

C  THE  RATE  COEFFICIENTS  FOR  THE  REVERSE  (EXCITATION)  REACTIONS  ARE  IN  THE 
C  FORM: 

C  RP(I)=01(I)*EXP(-F(I)/T)*R(I). 

C  HERE,  T  MAY  BE  EITHER  THE  NEUTRAL  OR  THE  ELECTRON  TEMPERATURE. 

DATA  A/1.68E-6,1.72E-5,1.52E-6,7.5E-10,6.2E-10,5.E-ll,l.E-9,l.E-9, 

I  l.E-10,l.E-9,l.£-10/ 

DATA  8/3*-. 5, 8*0./ 

DATA  C/37405.5, 105870. 9,267182.1, 8*0./ 

DATA  0/3*0. 2, 8*0./ 

DATA  8/3*0.356,8*1./ 

DATA  Dl/2.,l.,l.,3*2.,5*l./ 

DATA  F/12997. 0,36786. 3, 92836. 0,3*12997.0, 3*36786.3, 2*92836.0/ 

C  fV)(I)  ARE  EINSTEIN  'A*  COEFFICIENTS  (SEC-1)  (W)  RBI  1)  ARE  THE 
C  CORRESPONOIN6  RHO*B  FAaORS  (SEC-l)  FOR  SUNLIGHT  CONDITIONS. 

DATA  (VV6.£4,1.4E7,2.E7/ 

DATA  RB/8.5E-2,.17,2.1E-5/ 

DO  10  1=1,3 

R(  1 )  =A(  I  )*n**8(  I  )*(  D(  I  )+Q(  1  )*EXP(  -C(  1  )/TT) ) 
RP(I)=01(I)*EXP(-F(I)/TT)*R(I) 

10  CONTINUE 
DO  20  1=4,11 

R( I ) =A( I )*T**8( I )*( D( 1 )+G( I )*EXP( -C( I )/T) ) 
RP(I)=D1(I)*£XP(-F(I)/T)*R(I) 

20  CONTINUE 
DO  30  J=l,3 
I=J+11 
R(I)*^¥»(J) 

RP(I)=FLOAT(K)*RB(J) 
bO  CuriTIHuE 

C  0N2P(1)  IS  THE  DELTA  FAQOR  FOR  THE  GROUND  ELECTRONIC  STATE  doublet 
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C  SIGW  5+)  OF  N2+.  IT  IS  APPROXimTELY  EQUAL  TO  ONE-HALF  THE  ELECTRONIC 
C  PARTITION  FUNCTION. 

CN2P(l)=QEN2P(T)/2. 

C  the  delta  factors  CN2P(I)  (1=2  TO  4)  REFER,  RESPECTIVELY,  TO  DELTAS  FOR 
C  THE  3  STATES:  A  doublet  PI  u,  S  doublet  SI91A  u-t-,  ftit)  C  doublet  SIGMA  u-F. 
E=0. 

I F( K . EQ. 1 )E=£XP( 12997 .0/T+ALOG( .5*RB( 1) ) ) 

FAC=EXP(  12997 . 0*(  1  ./T-1  ./TT) ) 

DN2P( 2) =( R( 1 )*XNE*FAC+R( 4)*XN2+R( 5)*02+R( 6)*(XN+0)+E)/ 

1  (R(1)*XNE+R(4)*XN2+R(5)*02+R(6)*(XNHJ)+R(12))*0N2P(1) 

£=0. 

I F ( K . EQ . 1 ) E=£XP( 36784 . 0/T+AL06( RB( 2) ) ) 

FAC=€XP( 36784 . 0*( 1 ./T-1 ./TT) ) 
ON2P(3)=(R(2)*XNE*FAC+R(7)*XN2+R(8)*02+R(9)R(XN+0)+E)/ 

1  (R(2)*XN£+ft(7)*XN2+R(8)*02«(9)*(XNfO)+R(13))*ON2P(l) 

E=0. 

I F( K . EQ. 1 ) E=£XP( 92862 . 2/T+AL0G( RB( 3) ) ) 

FAC=£XP( 92862 , 2*( 1 ./T-1 ./TT) ) 

0N2P(  4) =(  R(  3)*XNE*FAC+R(  10)*()(N2+02)+R(  11  )lt()(NtO)+E)/ 

1  (R(3)*)(NE+R(10)*()(N2tO2)+R(ll)*()(N«))+R(14))*0N2P(l) 

RETURN 

END 


SUBROUTINE  OEL02()(NE,T,n,K,D02) 

C  THIS  SUBROUTINE  COMPUTES  DELTA  FACTORS  FOR  THE  VIBRATIOM^L  STATES  V'=0 
C  TO  20  OF  ELECTRONIC  STATE  02(8  TRIPLET  SIGMA  U  MINUS)  .THESE  DELTA 
C  FACTORS  CAN  BE  USED  TO  CONVERT  LTE  EMISSION  RATES  FOR  THE  SCHimW- 
C  RUN6E  aAN0S(8jV'->XiV")  TO  NLTE  RATES. 

DIMENSION  QP(24,21),A(21,24),F(21,24),Q(21,24),XNU(21,24),E(24), 

1  RB(21,24), 002(21), XK(21,24) 

DATA  QP/3.53E-9,l.fl2E-7,1.43E-6,1.28E-5,8.26E-5,4.09E-4,1.61E-3, 

1  5.15E-3,1.36E-2,3.04E-2,5.77E-2,9.28E-2,1.28E-1,.15,.152,.132, 

2  9.97E-2,6.6E-2,3.82£-2,1.93E-2,8.48E-3,3.24£-3,1.13E-3,3.81E-4, 

3  3.49E-8,9.09E-7,1.14E-5,9.07E-5,5.11E-4,2.17E-3,7.16E-3,1.88E-2, 

4  3.91E-2,6.54E-2,8.64€-2,8.59E-2,5.85E-2,1.98E-2,1.97E-5,1.95E-2, 

5  6.82E-2,.114,.131,.115,8.11E-2,4.73E-2,2.38E-2,1.09E-2,2.14£-7, 

6  5.06E-6,5.7E-5,4.E-4,1.96E-3,7.11E-3,1.95E-2,4.09E-2,6.46E-2, 

7  7.48E-2,5.77E-2,2.17E-2,9.87E-5,1.82E-2,5.76E-2,7.12E-2,4.13E-2, 

8  4.67E-3,8.84E-3,5.54E-2,.104,.119,.lfl2,7.29E-2,9.14E-7,1.96E-5, 

9  1.98E-4,1.23E-3,5.24E-3,1.61E-2,3.62E-2,5.9E-2,6.62E-2,4.49E-2, 

1  1.08E-2,1.7E-3,2.95E-2,5.38E-2,3.59E-2,3.27E-3,1.16E-2,5.29E-2, 

2  6.42E-2,2.84E-2,1.86E-4,2.27E-2,7.71E-2,.125,2.91E-6,5.68E-5, 

3  5.16E-4,2.84E-3,l.(J5E-2,2.71E-2,4.9E-2,5.94E-2,4,17E-2,9.66E-3, 

4  2.09E-3,2.86E-2,4.4E-2,1.96E-2,9.66E-5,2.5S-2,5.E-2,2.56E-2, 

5  1.65E-6,2.71E-2,6.18E-2,4.46E-2,5.58£-3,1.26E-2,7.66E-6,1.36E-4, 
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6  1.1IE-3,5.-I2E-3,1.’3£-2,3.71£-2,5.26£-2,4.37E-2,1.39E-2,4.5E-4, 

7  2.22E-2,3.85E-2,l.e5E-2.4.13E-4,2.57E-2,3.31E-2,1.04E-2,4.61E-3, 
3  3.36£*2,4.05E-2,5.22E-3,l.a8E-2,5.02E-2,5.06£-2,1.73£-5,2.81£-4, 
9  2.t37E-3,8.92E-3,2.45E-2,4.32E-2,4.6E-2,2.26E-2,5.67E-4.1.26E-2, 

1  3.33E-2.1.39E-2.1.1E-5,1.33E-2,3.31£-2,7.8£-3,5.86£-3,3.55E-2, 

2  2.47E-2,2.E-b,2.61E-2,4.44E-2,1.29E-2,4.86E-3,3.5E-5,5.2E-4, 

3  3.46E-3,1.32E-2,3.09E-2,4.42£-2,3.33E-2,6.55E-3,3.41£-3,2.6E-2, 

4  2.55£-2,2.22E-3,1.02E-2,2.95£-2,1.08£-2,2.83E-3,2.88£-2,2.05£-2, 

5  9.68£-5,2.6£-2,3.26E-2,2.57£-3,1.49£-2,4.15E-2,6.3£-5,8.59£-4, 

6  5.17E-3,1.74E-2,3.49E-2,3.94E-2,1.88E-2,1.4£-4,1.36E-2,2.74E-2, 

7  9.54E-3,2.03£-3,2.29E-2,1.68£-2,5.78E-5,2.01£-2,2.18£-2,1.35E-4, 

8  1.91£-2,2.7E-2,1.27E-3,1.69£-2,3.36£-2,5.26E-3,1.01£-4,1.27E-3, 

9  6.94E-3,2.07£-2,3.53E-2,3.fl7E-2,7,57£-3,2.fl9E-3,2.09E-2,1.91E-2, 

1  6.35£-4,1.15E-2,2.13E-2,2.66£-3,8.8E-3,2.28£-2,3.44£-3,9.29E-3, 

2  2.48E-2,3.39E-3,l.llE-2,2.88£-2,5.04E-3,l.£-2,1.48£-4,1.72E-3, 

3  8.61£-3,2.29E-2,3.31£-2,2,15E-2,1.59E-3,7.53E-3,2.2E-2,8.8E-3, 

4  l,39E-3,1.79£-2,1.13E-2,6.04£-4,1.76£-2,1.12E-2,1.16E-3,2.£-2, 

5  9.57E-3,3.07E-3,2,39E-2,8.44£-3,4.57£-3,2.54E-2,1.96E-4,2.14£-3, 

6  9.79E-3,2.33£-2,2.85E-2,1.31£-2,1.06E-7,1.22E-2,1.78E-2,2.1E-3, 

7  6.44E-3,l.E9E-2,2.77E-3,6.12E-3,1.65E-2,1.6l£-3,9.07E-3,1.61E-2, 

8  1.84E-4,1.43E-2,1.43E-2,2.64£-4,1.9E-2,1.17E-2,2.4£-4,2.45£-3, 

9  1.04E-2,2.23£-2,2.31£-2,7.13E-3,9,34£-4,1.4X-2,1.19E-2,2.05£-5, 

1  1 .05£-2,l 46£-2,l .36E-5,1 .09E-2,9.79E-3,3.26E-4,1.38E-2,6.59E-3, 

2  2.86E-3,1.64E-2,2.55E-3,8.23E-3,1.72E-2,3.83£-4,2.7E-4,2.6£-3, 

3  1.03E-2,2.£-2,1.77E-2,3.36E-3,2.57E-3,1.36E-2,6.77E-3,6.81E-4, 

4  1.15E-2,6.19E-3,1.06E-3,1.16E-2,3.78E-3,3.27E-3,1.22E-2,9.51E-4, 

5  8.04£-3,1.04E-2,1.83E-4,1.37E-2,6.92E-3,2.69E-3,2.81£-4,2.58E‘3, 

6  9.58E-3,1.17£-2,i.3E-2,1.34E-3,3.81£-3,1.16E-2,3.33E-3,2.05E-3, 

7  1 .03£-2,2.57E-3,2.95E-3,9.55E-3,7.7E-4,5.83£-3,7.99E-3,8.04E-5, 

8  9.81£-3,4.18E-3,3.02E-3,1.2E-2,8.9E-4,7.74£-3,2.72E-4,2.4E-3, 

9  8.47E-3,1.4E*2,9.29E-3,4.29E-4,4.33E-3,9.06E'3,1.43£-3,3.0l£-3, 

1  3.08E-3,7.94E-4,4.11E-3,6.74£-3,7.23E-6,6.55E-3,4.25£-3,1.16E-3, 

2  8.5iE-3,9.S9E-4,5.45E-3,7.48E-3,1.15E-4,9.11E-3,2,48E-4,2.12E-3, 

3  7,15E-3,1.12E-2,6.52E-3,9.27E-5,4.21E-3,6.75E-3,5.29E-4,3.32E-3, 

4  5,9E-3,1.42E-4,4.32E-3,4.34E-3,1.81E-4,5.92E-3,1-3«-3,2.1*-3, 

5  6,18E-3,4.04E-5,5.98E-3,3.76E-3,1.23E-3,7.61£-3,2.13E-4,1.77E-3, 

6  5.76E-3,8.52E-3,4.51E-3,5.82E-6,3.69E-3,4.84E-3,l.ElE-^3.11E-3, 

7  4,l£'3,1.67E-6,3.87E-3,2.6SE-3,4.93E-4,4.72E-3,7-72E-4,2.56E-3, 

8  4.06£-3,8.12E-5,5.22E-3,l.«E-3,2.13E-3,5.32E-3,1.74E-4,1.41£-3, 

9  4.5E-3,6.«-3,3.1£-3,3.49E-6,3.03E-3,3.41E-3,3.63E-5,2.64E-3, 

1  2.79E-3,2.4£-5,3.17E-3,1.59E-3,6.64E-4,3.52E'3,2'74E-4,2.41£-3, 

2  2.55E-3,2.94E-4,4.06E-3,6.21E-4,2.73E-3,3.43E-3,1.35£-4,1.08E-3, 

3  3.37E-3.4.6«-3,2.11E'3,1.69E-5,2.34E-3,2.34E-3,3.88E-6,2.09E-3, 

4  1.87E-3,6.36£-5,2.42E-3,9.46£-4,6.74E-4,2.K-3,8.45£-5,2.01£-3, 

5  1.57E-3,4.14£-4,2.94E-3,2.13E-4,2.14£-3,2.11£-3,8.81£-5,6.99E-4, 

6  2.15E-3,2.91E-3,1.25E-3,2.13E-5,1.52E-3,1,4E-3,3.43E-8,1.37E-3, 

7  1.09E-3,7.E-5,1.56E-3,5.03E-4,5.11£-4,1.52E-3,2.05E-5,1.35£-3, 
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3  3 . 48E-4 , 3 . 64E'4 , 1 . 8E-3 , 6 . 05E-5 , 1 . 51E-3 , 1 . 14E-3/ 

C  >;!BRATirm  PARPHETERS  FOR  02(X) 

WTA  M.4)(,WY,WZ/1580. 193, 11. 381,. 04747, -.001273/ 

DO  20  1=1,21 
OP=FL»T(I-l) 

DO  10  J=l,24 
VPP=F'wOAT(J-l) 

C  WWENUHeERS  FOR  THE  SDOm^-RUNGE  8PNOS  (V'-;v".'' 

XNU(I,J)=49793.29f709.31*(VP+.5)-10.65*(VW-.5)**2-.139*(W.5)**3 
1  -(W*(OPW'.5)-W<*(yPP+.5)**2+WY*(WH.5)**34WZ*(VPP+.5)**4) 

C  FIVWCK-CONOON  FACTORS  FOR  THE  SCHIJMAM4-  RUN6E  BANOS: 

Q(1,J)=QP(J,I) 

C  ABSORPTION  OSCILLATOR  STRENGTHS, F (I, J),0£FMD  IN  TERMS  OF  THE  FRANCK- 
C  CONOCN  FACTORS: 

F(I,J)=2.36£-6*Q(I,J)*XNU(1,J) 

C  EINSTEIN  A  COEFFICIENTS  (SEC-1)  DEFINED: 

A( I , J)=0 .667*XNU( I , J)**2*F( I , J) 

C  RHO*B  FACTORS  (SEC-1)  FOR  SIRIGHT  EXCITATION  OF  02(8;y')  DEFINED: 

R8( I , J) =5 . 4E-6*A( I , J)/( EXP( 2 , 391E-4*XNU( 1 , J) ) -1 . )*FL0AT( R) 

C  RATE  COEFFICIENT  FOR  ELECTRON  EXCITATION  OF  THE  STATE  02(B;V')  FROM 
C  THE  STATE  02(Xiy*): 

Xl((  I ,  J)=.754*F(  I ,  J)/XNO(  I ,  J)**l  .5*EXP(  .3366*(n/)WJ( ! .  J)  )«.546- 
1  1.4388*)(NU(I,J)/7r) 

10  CONTINUE 
20  CONTINUE 

C  DC  IS  THE  GROUND-STATE  DELTA  FACTOR 

00=1  .+1 . 336E-2*T** .  445«XP( -Ifl309.1/T) 

00  40  1=1,21 
S1=0. 

S2=0. 

S3=fl. 

DO  30  J=l,24 

S1=S1+XK( I , J)*EXP(1,4388*XNU( I , J)/TT) 

S2=S24A(I,J) 

S3=S3+(  RB(  1 ,  J)+)^*XK(  I ,  J)  )*£XP(1 . 4388»XNU(  1 ,  J)/T) 

30  C0NT1N(£ 

FN=S3*00 

F0=S2tXNE*Sl 

C  DELTA  FAaORS  FOR  THE  STATES  02(B;V'=0  TO  20)  DEFINED: 

002(I)=FN/FD 
40  CONTINUE 
RETURN 
END 
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SUBRaiTINE  DEL>«(XN2,02,XNO,:02,H20,M,T,n,{)N2.K,CNO) 

C  ’■HIS  SUBPOlfTINE  COMPUTES  DELTA  FACTORS  FOR  7  ELECTRONIC  STATES  OF  NO. 
DIMENSION  R(27'i,RP(27),f^(ll),RB(ll),DNO(7),A(b),S(6),C(6),D(6), 

1  XKtlO),DN2(ll) 

C  ELECTRONIC  PARTITION  FLNCTIQN  FOR  N2. 

0EN2(X)=l.+3.*O(P( -71587. 2A)+6.*(EXP(-B5330.6A)+£XP(-85435.9A)) 

1  +3. *EXP( -94748. 4A)+EXP( -97463. 3A)+2.*(EXP(-99207.0A)+EXP(-1031 

2  59. BA) )+5.*EXP(-112082.5A)+6.*(EXP(-125319.5/'X)+EXP( -128021. 3A 

3  ))+10.*EXP( -134096. 2A)+3.*EXP(-137800.4A)+6.*EXP( -140432. 0A)+ 

4  EXP( -142211. 4A)+2.*EXP(-145056.2A) 

C  ELECTRONIC  PARTITION  FINCTION  FOR  NO. 

OI£NO(X)=4.*(1.+.5*£XP( -65119. 5A)+£XP(-71651.6A)+EXP( -75349. 3A)+ 
1  .5*EXP( -81291. 5A)+EXP(-91938.5A)+.5*EXP( -92657. 9A)) 

C  TWO-BOOY  RATE  COEFFICIENTS  FOR  ELECTRON-IMPACT  OEEXCITATION  ARE  IN 
C  THE  FORM!  R(I)=A(I)*£XP(B(I)*n**.546).  THE  COSUIESPCNDING  ELECTRON- 
C  IMPACT  EXCITATION  RATE  COEFFICIENTS  ARE  IN  THE  FORM:  RP(I)<(I)*EXP( 

C  -0(I)/n)*R(I)  WHERE  n  IS  THE  ELECTRON  TENTERATUREO  DEGREES  KELVIN). 
DATA  A/9 . 6E-11 , 9 . 5E-13 , 3 . 5E-10 , 4 . 8E-10 , 5 . lE-11 , 1 . 7E-10/ 

DATA  B/9.66£-4,9.17E-4,8.93E-4,8.56£-4,8.E-4,7.97E-4/ 

DATA  C/.5,l.,l-.-5,l.,.5/ 

DATA  0/65119. ,71652. ,75349. ,81291 . ,91938. ,92658./ 

C  (¥)(!)  ARE  THE  EINSTEIN  ‘A*  COEFFICIENTS  (SEC-1)  AND  RB( I)  ARE  THE 
C  CORRESPONDING  RHO*e  FACTORS  (SEC-1)  FOR  SUNLIGHT  CONDITIONS. 

DATA  (VV'1.7E6,2.3E4,1.E7,1.8E7,3.£6,1.E7,1.E7,1.3E7,9.3E5,2.9E5, 

1  1.6E7/ 

DATA  RE/9.2£-5,8.4E-7,2,£-4,6.6£-5,3.8E-6,5.5E-6,18.,8.9,9.4E-2, 

1  .12,17./ 

C  .XK  ARE  rate  COEFFICIENTS  FOR  COILISIONAL  QUENCHING  BY  HEAVY  PARTICLES 
DATA  XK/1.4E-13,6.9E-11,3.E-12,7.2E-12,3.7E-11,6.E-10,5.E-10, 

1  4.5C-11,1.4E-11,5.E-11/ 

DO  10  1=1,6 

R(I)=A(I)*EXP(B(I)*n«.546) 

RP(I)=C(I)*EXP(-0(I)/TT)*R(I) 

10  CONTINUE 
DO  15  1=7,16 
J=l-6 
R(I)=XR(J) 

15  CONTINUE 

RP(ll)=.3«XP(-3763./T)*R(ll) 

C  OTHER  COLLISIONS.  QUENCHING  RATE  COEFFICIENTS. 

R 702=1, 5E-10 

R)NO=2.E-10 

R7C02=3.8E-10 

R.7H20=7.7E-lfl 

R1002=5.4E-11 

R10NO=2.£-10 

R10H20=1.3E-10 

R14O^4.3E-10 
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R14NO=5.E-10 
R14C02=6.3£-10 
DO  20  1=17,27 
J=l-16 
RU  )=*¥>(  J) 

RP(I)=FLOAT(Kj*ftB(J) 

20  CONTINUE 

C  CNO(l)  IS  THE  DELTA  FACTOR  FOR  THE  GROUND  ELEaRONIC  STATE  (X  doublet 
C  PI)  OF  NO.  IT  IS  APPROXIMATELY  EQUAL  TO  ONE-QUARTER  OF  THE  ELECTRONIC 
C  PARTITION  FUNCTION. 

ONO(l)=Q£NO(T)/4. 

81=R(l)*XNE+R(7)*XN2+R702*02+R7NO*XNOtR7C02*C02+R7H20*H20FR(17)+ 

1  RP(23)+RP(24)+RP(25) 

82=2. *(R(12)*XN2+R( 23) )*EXP( -10230 ./T) 

83=EXP(AL06(R(24))-16172./T) 

84=EXP(AL06( R( 25) ) -27538 ./T) 

E=0. 

1F( K . EQ . 1 ) E=£XP( 65113 ./7+AL0G( RB( 1 )) ) 

85=2.*(.5*£XP(65119.*(1./T-1./TT))*R(1)*XNE+E+3.*R(8)*XN2*0N2(2)* 

1  EXP(-6467./T)/Q£N2(T)) 

B6=RP(23) 

87=2.*(R(  3)*XNE+R(11  )*XN2*0N2(  1  )/Q£N2(T)+R(12)*»C+R(  13)*(02+XN0) 

1  +R( 19)+R( 23)+RP( 26) )*EXP( -11741 ./T) 

B8=EXP(AL06(R(26) ) -27538. /T) 

E=0. 

IF(K.EQ.l)£=O(P(65119./r+AL0G(RB(3))) 

B9=2.*(  EXP(  65119./T-75349./TT)*R(  3)*XNE+E+3.*RP(  11  )*XN2*l)N2(  2)* 

1  EXP(-6467./T)/QEN2(T)) 

B10=RP(24) 

B11=£XP(AL06(R(4)*)(NE+R(14)*XN2+R1402*02+R14NO*XN(HR14C02*C02+ 

1  R(20)+R(24)+RP(27))-16172./T) 

B12=£XP(AL06(R(27) )-27538./T) 

E=0. 

1 F ( K . EQ . 1 ) E»£XP( 6511 9 ./T+AL06( RB( 4) ) ) 

813=2.*( .  5*€XP(  65119./T-81291  ./TT)*R(4)*XNE+E) 

B14=RP(25) 

815=0. 

IF(K.£0.1)815»€XP(ALOG(2.*R8(10))-10230./r) 

816=0. 

I F  ( X .  EQ .  1 )  B16=£)(P(  ALOG(  RB(  1 1 ) ) -1 6172  ./T) 
B17=EXP(AL0G(R(6)*XNE+R(16)*XN?fR(22)+R(26)+R(27))-27538./T) 

E=0. 

I F ( K . EQ . 1 ) E=£XP( 651 13 ./T+AL06( R8( 6) ) ) 

81 8=2 ,  R( .  5*EXP(  6511 9  ./T-92658  ./TT)*R(  6)*)(NE+E) 

G=81 0*( B3*( 88*815-87*81 7) -816*1 84*87+82*88) ) -Bll*( Bl*( 88*815- 

1  B7*817)+B2*(B8*814+86*817)+84*(B6*8l5tB7*8l4))-B12*(87*(Bl*816+ 

2  B3*814)+B6*(B3*B15-B2*816)) 

A1=85*CN0(1) 
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A2=89tW0(l) 

A3=813*CN0(1) 

A4=818*ONO(1) 

C  THE  DELTA  FACTORS  i>JO(I)  il=2  TO  7)  REFER,  RESPECTIVELT,  TO  DELTAS  FOR 
C  THE  STATES:  A  doublet  SIGMA  +,  8  doublet  PI,  C  doublet  PI,  D  doublet 
C  SIGMA  +,  8'  doublet  DELTA,  AM)  E  doublet  SIGMA  +. 

DNO( 2) =t B2*( A2*( Bll*817-B12*816)+88*( A3*81&+A4*811 ) )+B15*(  B3*( A2* 

1  812-A3*B8)+B11*(A2*B4-A1*88))+B7*(A1*(B11*817-B12*B16)4A3*( 84*816 

2  +B3*617)+A4*(83*ei2+B4*B11)))/G 

DNO(3)=((RP(2)*XNE+RP(18))*EXP(71652./T)+3.*R(9)*XN2*ON2(2)/Q£N2(T 

1  ))*DM)(l)/(R(2)*XNE+R(10)*){N2+R1002*02+R10NO*XNO+R10H20*H20tR(18> 

2  ) 

DNO(  4)=(  B6*(A1*(  811*81 7-B12*B16)+A3*(B3*B17+64*616)-HA4*(  B3*812+B4* 

1  811))+A2*(81*(811*fll7-B12*816)-B10*(83*817+B4*816)-S14*(B3*612+B4 

2  *B11))+B8*(B1*(A3*B16+A4*B11)+B10*(A1*616-A4*63)+B14*(A1*B11+A3* 

3  B3)))/6 

0NO(5)=(810*(B2*(A2*B17+A4*88)+B7*(Al*8174A4*B4)+B15*(A2*B4-Al*B8) 

1  )-A3*(Bl*(B8*BlS-B?*817)+B6*(B2*817+84*815)+B14*(B2*B8+fi4*B7) )+ 

2  B12*(  Bl*( A4*B74A2*615)  -B6*( A4*82-Al*B15)+614*(A2*82-hAl*87) )  )/S 
DM)(  6)=(  RP(5)*XNE+RP(  21 )  )*EXP(  91938. /T)*ONO(1 )/( R(  5)*)(Ne+R(  15)*XN2 

1  W21)) 

ONO( 7) s( B6*( B15*( A3*B3+Al*eil ) -e2*( A4*eil+A3*B16) )+67*( Bl*( A4*811+ 

1  A3*B16)  -B10*(  A4*83-A1*B16)+014*(  A3*B3W1*811 )  )+A2*(  B2*(  B10*B16+ 

2  811*814)-B15*(B3*B10-Bl*«l)))/6 
RETURN 

END 


SUBROUTINE  0eLNOP()tfC,T,TT,K,O0,01,0NOP) 

C  THIS  SUBROUTINE  COMPUTES  DELTAFAGORS  FOR  THE  VIBRATIONAL  STATES  v'=0 
C  TO  10  OF  ELECTRONIC  STATE  NOf(A  singlet  PI),  THESE  DELTA  FACTORS  CAN 
C  BE  USED  TO  CONVERT  LTE  EMISSION  RATES  FOR  THE  BAER-MIESCHER  BANOS 
C  (A;v'-->X;v‘)  TO  NLTE  RATES. 

DIMENSION  g(ll,8),A(ll,8),XNU(ll,8),RB(ll,8),XK(ll,8), 

1  ONOP(ll) 

C  ELECTRONIC  fWTITIQN  FUNCTION  FOR  NOF 

0£NOP(X)=1.+3.*(E)(P(-75090.A)+2.*E)(P( -85234  ,A)+2.*E)(P(-89032./X) 

1  tC<P(-97435.A)) 

C  TX  TA  («£  TIC  AflSaUTE  ENERGIES  (CM-1)  OF  VIBRATIONAL  STATES  OF 
C  THE  ELECTRONIC  STATES  (X  singlet  SIOIAf)  fiNO  (A  singlet  PI),RESPECTIVaY. 
rX(X)=2376.42*(X+.5)-16.262*(X+.5)**2-.01562*(X+.5)**3 
TA(X)=73083. 5+1601, 93*(X+.5)-20.207*(X+.5)**2-.2683*(X+.5)**3 
DATA  Q/,D217,. 0643,, 1068,, 1321,, 1363,, 1246,, 1047,, 083,, 063,, 0465, 

1  ,0336, ,0962, ,1617, .1338, .0644, ,0138, 9.024E-5,. 0121, .0323, .049, 

2  .0581,. 0596,. 1958,. 1329,. 0158,. 2101, .0555,. 0744,. 0569,. 0274, 

3  .0642,3.1E-7,4.598E-3,.2443,.0214,.0369,.0922,.0507,4.54£-3, 


.\-ll 


4  b.:33E-3,3.258£-2,4.857E-2.4.525E-2,3.037E-2,.21,.fll81,.113, 

5  . 02764, 7. 03E-3,. 0532,. 0577,. 0247, 1.66E-3, 4. 31E-3,. 0198,. 133, 

6  .1205, .05, .0187, .078, .0327, 4. 31E-5, .024, .0462, .0378, .0157, .0644, 

7  . 1835, 1.47E-3,. 0963,. 0214,. 0141,. 0579,. 0371, 3. 32E-3, 5. 96E-3, 

8  .0266, .0246. .1549, .0788, .0509, .0209, .0712, .0142,9.9196*3, .0419, 

9  .0392,. 0123/ 

DO  20  1=1,11 
'7P=FLaAT(I-l) 

DO  10  J=l,8 
VPP=EL0AT(J-1) 

C  WAVENUMBERS  FOR  THE  BAER-MIESCHER  BANOS  (v'->v')  : 

/NU(I,J)=TA(VP)-'n((VPP) 

C  EINSTEIN  A  COEFFICIENTS  (SEC-l)  DEFINED  : 

A( I , J)=6.429E-8*XNU(I , J)**3*0( I , J) 

C  RHO*B  FAaORS  (SEC-1)  FOR  SUMKiHT  EXCITATION  OF  NO+(A;v')  DEFINED  : 

m i , J)=l .08£-5*A( I , J)/(EXP(2.391£-4*XNU( I , J) )-l .)*fL0AT(K) 

C  RATE  COEFFICIENT  FOR  ELECTRON  DEEXCITATION  OF  THE  STATE  NO+(A:v')  TO 
C  THE  STATE  NOf(X;v*)  : 

XR( I  ,J)=20.44*A(  I ,  J)/(SQRT(n)*XNU(l  ,J)**3)*(  .2+.356*EXP(-4.14* 

1  XNU(I,J)/TT)) 

10  CONTINUE 
20  CONTINUE 
DO  50  1=1,11 
S1=0. 

S2*0. 

DO  30  Js2,8 

£=0. 

IF(K.EQ.l)£=£XP(ALCI6(RB(I,J)/2.)+l.  4388*XNU(  1 ,  J  )/T ) 

Sl=Sl+(  E+)(NE*EXP(  1 . 4388*XNU(  I ,  J)*(  1  ./T-1  ./TT)  )*XK(  I ,  J) ) 

1  *(01/00 )**(J-2) 

30  CONTINUE 
DO  40  J=l,8 

S2=S2+A(I,J)+XNE*XK(I,J) 

40  CONTINUE 

£=0. 

iF(K,EQ.l)E=EXP(AL06(R8(I,l)/2  )+1.4388*XNU(I ,1)/T) 

FN=:(  E+XNE*EXP(  1 . 4388*XNU(  1 ,1  )*(1  ./T-1  ./TT)  )*XK(  1,1)  )*00+01*S1 
ONOP(I)=FmO£NOP(T)/S2 
50  CONTINUE 
RETURN 
ENO 
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